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Abstract 

Motivated by physical and numerical observations of time oscilla- 
tory "galloping" , "spinning" , and "cellular" instabilities of detonation 
waves, we study Poincare-Hopf bifurcation of traveling- wave solutions 
of viscous conservation laws. The main difficulty is the absence of a 
spectral gap between oscillatory modes and essential spectrum, pre- 
venting standard reduction to a finite-dimensional center manifold. 
We overcome this by direct Lyapunov-Schmidt reduction, using de- 
tailed pointwise bounds on the linearized solution operator to carry 
out a nonstandard implicit function construction in the absence of a 
spectral gap. The key computation is a space-time stability estimate 
on the transverse linearized solution operator reminiscent of Duhamel 
estimates carried out on the full solution operator in the study of 
nonlinear stability of spectrally stable traveling waves. 
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1 Introduction 

Motivated by physical and numerical observations of time-oscillatory "gal- 
loping" or "pulsating" instabilities of detonation waves |MTt IBMR} IFW| 
IMTt I ALT} IAT| \F1\ IF2] , we study the Poincare-Hopf bifurcation of viscous 
shock waves in one spatial dimension. Our main result is to obtain a rigor- 
ous (nonlinear) characterization in terms of spectral information. The com- 
plementary problem of verifying this spectral information has been studied 
already in |LyZl , LyZ2| , where it was shown that transition from stability 



to instability of viscous detonation waves generically involves a Poincare- 
Hopf type bifurcation in the spectral configuration of the linearized operator 
about the wave. By essentially the same analysis we obtain also a corre- 
sponding multi-dimensional result applying to planar viscous shock fronts 
traveling in a cylinder of finite cross-section, with artificial Neumann or pe- 
riodic boundary conditions. This gives a simplified mathematical model for 
time-oscillatory instabilities observed in detonation waves moving within a 
duct, which, besides the longitudinal galloping instabilities described above, 
include also transverse "cellular" or "spinning" instabilities. The method 
of analysis appears to be of general application, in particular, with suitable 
elaboration, to extend to the originally motivating case of viscous detonation 
waves of the reactive Navier-Stokes equations with physical viscosity. 

From a mathematical standpoint, the main issue in our analysis is that 
the linearized equations about a standing shock wave have no spectral gap 
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between convective modes corresponding to essential spectra of the linearized 
operator L about the wave and oscillatory modes corresponding to pure imag- 
inary point spectra. This prohibits the usual PDE analysis by center man- 
ifold reduction to a finite-dimensional subspace; likewise, at the linearized 
level, decay to the center subspace is at time-algebraic rather than time- 
exponential rate, so that oscillatory and other modes are strongly coupled. 
We overcome this difficulty by the introduction of a nonstandard Implicit 
Function Theorem framework suitable for infinite-dimensional Lyapunov- 
Schmidt reduction in the absence of spectral gap (Section 12.1 1) , augmented 
by a finite-dimensional "weak" Implicit Function Theorem suitable for finite- 
dimensional bifurcation analysis in situations of limited regularity (Section 
12.31) . The latter, based on the Brouwer Fixed-point Theorem rather than 
the standard Contraction-mapping construction, seems of interest in its own 
right. 

The study of bifurcation from stability is a natural followup to our pre- 
vious work on stability of viscous shock and detonation waves; see, e.g., 
[ZHl IMaZ3l IMaZ4l iHRZl [LyZTj [LyZ2l ILRTZj . Interestingly, the key estimate 



needed to apply our bifurcation framework turns out to be a space-time sta- 
bility estimate on the transverse linearized solution operator (meaning the 
part complementary to oscillatory modes) quite similar to estimates on the 
full solution operator arising in the stability analysis of viscous traveling 
waves. To carry out this estimate requires rather detailed pointwise infor- 
mation on the Green function of the linearized operator about the wave; see 
Sections ll.6.4l andl5l Indeed, we use the full power of the pointwise semigroup 
techniques developed in [ZHt IMaZl[ IMaZ3j . 



1.1 Shocks, detonations, and galloping 

The equations of compressible gas dynamics in one spatial dimension, like 
many equations in continuum mechanics, take the form of hyperbolic con- 
servation laws 

(1-1) u t + f(u) x = 

(Euler equations) or hyperbolic-parabolic conservation laws 

(1.2) u t + f(u) x = {B{u)u x ) x 

(Navier-Stokes equations), depending whether second-order transport effects- 
in this case, viscosity and heat-conduction- are neglected or included. Here, 



4 



x, t G M 1 are spatial location and time, u G W 1 is a vector of densities of 
conserved quantities- mass, momentum, and energy in the case of gas dy- 
namics, / G M. n is a vector of corresponding fluxes, and B G W ixn is a matrix 
of transport coefficients. 

Such equations are well known to support traveling-wave solutions of form 



'1.3) u(x,t) = u(x — st) :- 



(discontinuous) and 



w_ x — st < 0, 
u + x — st > 



(1.4) u(x, t) — u(x — st) , lim u(z) = u± 

(smooth), respectively, known as ideal and viscous shock waves. These waves 
may be observed physically; indeed, the corresponding physical objects ap- 
pear to be quite stable [BEJ . 

Similarly, the equations of reacting, compressible gas dynamics take the 
form 

^ u t + f(u, z) x - qK(j){u)z = 

z t + (v(u, z)z) x + Kcj){u)z = 

(reactive Euler, or Zeldovich-von Neumann-Doering equations (ZND)) or 

u t + f(u, z) x - qK(j)(u)z = (B(u, z)u x ) x , 

(1.6) 

z t + (v(u, z)z) x + K(p(u)z = (D(u, z)z x ) x 

(reactive Navier-Stokes equations (rNS)), where x, t G R 1 , u G 1R™ is as be- 
fore, z G W is a vector of mass fractions of different reactant species, v G M 1 
is fluid velocity, and /, B and D are flux vectors and transport matrices 
depending (through z) on chemical makeup of the gas, with f(u,0) = f(u), 
B(u,0) = B(u). The matrix K models reaction dynamics, q is a constant 
heat release coefficient (q > corresponding to an exothermic reaction), and 
<p G M 1 is an ignition function serving to "turn on" the reaction: zero for tem- 
peratures below a certain critical temperature and positive for temperatures 
above. 

Equations (jl.5p and (jl.6p support traveling wave solutions 
(1.7) u(x, t) — u(x — st) 
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analogous to (jl.3p and (11. 4p . known as ideal and viscous detonation waves, 
which may likewise be observed physically. However, in contrast to the shock 
wave case, detonations appear to be rather unstable. What is typically ob- 
served is not the planar, steadily progressing solution (11.71) . but rather a 
nearby solution varying time-periodically about (|1.7p : that is, an apparent 
bifurcation with exchange of stability. In the typical experimental setting of 
a detonation moving along a duct, these may be longitudinal "galloping", 
or "pulsating", instabilities for which the planar structure is maintained, 
but the form of the profile changes time-periodically, or they may be trans- 
verse instabilities in which the front progresses steadily in the longitudinal 
(i.e., axial) direction, but develops time-oscillatory structure in the trans- 
verse (cross-sectional) directions of a form depending on the geometry of the 
cross-section: "cellular" instabilities for a polygonal (e.g. rectangular) cross- 
section; "spinning" instabilities for a circular cross-section, with one or more 
"hot spots", or "combustion heads", moving spirally along the duct. 

Stability of shocks and detonations may be studied within a unified 
mathematical framework; see, e.g., (ErTJ ElS IKoTl IKo2l [0 [BTj EH ESI 
HI E [MD M [MSI EMI [Mil IGJLW] in the inviscid case flUO, flESD, and 
p iGoTllGo2lllMlKMmiMm iraP ICT 

IKK1 IZTj IZ2 j IZ31 IZlf IMaZ2l IMaZ3l [MaZil IGMWZll IGWMZ21 IGMWZ31 iHZl 
iBLl |LyZTl |LyZ2l iJLWl IHuZTl IHuZ2l IBHRZl IFZl IF31 ILRTZj in the viscous 
case (11.21) . (jl.fip . In particular, stability in each case has been shown to 
reduce to spectral considerations accessible by standard normal-modes anal- 
ysis (see Section [L2|) . However, they are studied with different motivations, 
and yield different results. Historically, it appears that instability motivated 
the physical study, which has focused on the detonation case and spectral 
stability criteria |BEj . By contrast, the mathematical study has focused on 
stability and the somewhat simpler shock wave case, with the main (diffi- 
cult!) issue being to establish full nonlinear stability assuming that spectral 
stability holds in the form of a suitable "Lopatinski" or "Evans" condition 
[M3llMel lZH lGWMZ2] . 

Regarding bifurcation, there is strong evidence that the detonation insta- 
bilities described above correspond to a Poincare-Hopf bifurcation (indeed, 
our description above benefits much by hindsight). In particular, Erpenbeck 
|Er2j and Bouiiicux. Majda, and Roytburd [BMRj have carried out formal 
asymptotics in support of this viewpoint for the one-dimensional (gallop- 
ing, or pulsating) case, in the context of the ZND equations. More recently, 
Kasimov and Stewart [KSJ have carried out a definitive study in the multi- 
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dimensional case by numerical linearized normal modes analysis, also for the 
ZND equations, in which they demonstrate that the onset of such instabili- 
ties indeed corresponds to crossing of the imaginary axis by a conjugate pair 
of eigenvalues toward the unstable (positive real part) side, with, moreover, 
excellent correspondence between observed nonlinear oscillations and the as- 
sociated normal modes. Depending whether the maximally unstable mode 
is in the longitudinal or transverse direction, the oscillatory behavior is seen 
to be of galloping or spinning (cellular) type. However, up to now, there has 
been carried out no rigorous analysis of these phenomena. To fill this gap is 
the object of the present study. 

1.2 Spectral considerations 

The main difficulty in the analysis of shock or detonation waves is the lack 
of a spectral gap between stationary or oscillatory modes and essential spec- 
trum of the linearized operator about the wave. In the inviscid, hyperbolic 
case (II. ip . (jl.Sp . the linearized problem has essential spectrum filling the 
imaginary axis. In the viscous case (ll.2p . (11.61) . the essential spectrum is, 
rather, tangent at the origin to the imaginary axis. 

For example, without loss of generality taking a standing-wave solution 
u = u(x), s = (i.e., working in coordinates moving with the wave), and 
linearizing (jl.2p about u, we obtain, taking B = I for simplicity, 

u t = Lu := u xx - (Au) x , A := df(u(x)), 

where the asymptotically constant-coefficient operator L is what we have 
called the linearized operator about the wave. A standard result of Henry 
( [He] . Theorem A. 2, chapter 5) on spectrum of operators with asymptoti- 
cally constant coefficients asserts that the rightmost (i.e. largest real part) 
envelope of the essential spectrum with respect to any LP is the envelope of 
the union of the rightmost envelopes of the spectra of the limiting, constant- 
coefficient operators at ±oo: in this case, 

(1.8) L ± := d 2 x - A ± d x , A±:=df(u±). 

The L? spectra of L± may be computed by Fourier transform to be the curves 
traced out by dispersion relations 

(1.9) A(0 := ia±£ - £ 2 , (Gt 1 , 
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where are the eigenvalues of A±. These eigenvalues are real and nonzero 
under the physical assumptions that (II. ip be hyperbolic in a neighborhood 
of u± and the shock be noncharacteristic. Likewise [S], L has always an L 2 
zero-eigenvalue with eigenfunction u'(x), associated with translation of the 
wave. 

In the absence of a spectral gap, standard stability and bifurcation theo- 
rems do not apply, and so we must carry out a refined analysis. 

1.2.1 Spatio-temporal description 

The spectral configuration just described translates in the x-t plane, in the 
stable case that there exist no other eigenvalues A of L in the nonnega- 
tive complex half-plane 3?A > 0, to the following description developed in 
|L1 | IL3 | [ZH| IZ4l IMaZ3j of the Green function G(x, t; y) associated with d t — L. 
A point source, or delta-function initial datum, originating at y > will prop- 
agate initially as an approximate superposition of Gaussians with constant 
mass (total integral), centered along hyperbolic characteristics dx/dt = 
determined by the asymptotic system at +oo. Those propagating in the 
positive direction will continue out to +oo; those propagating in the nega- 
tive direction will continue until they strike the shock layer at approximately 
x = 0, whereupon they will be transmitted and reflected along outgoing 
characteristic directions dx/dt = a^, < and dx/dt = a^, a\ > 0. 
In addition, there will be deposited at the shock layer a certain amount of 
mass in the stationary eigenmode u\x) (integral f u(x)dx = u + — w_ ^ 0), 
corresponding to translation of the background wave. 

We refer the reader to |ZHt IMaZ3j for further discussion and details. 
For the present purpose, it suffices to note that the Green function may be 
modeled qualitatively as the sum of terms 

(1.10) K(x, t; y) := r ^ e -^~v-^?l^ 
and 

(1.11) J{x, t; y) := u'(x)errfn((-y - at)/2t 1/2 ) 
propagating with noncharacteristic speed a < 0, where 

(1.12) errfn(z) := — 

2tt 
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the first modeling moving heat kernels (Gaussian signals), and the second 
the excitation of the zero-eigenf unction by incoming signals from y > 0. 

From this description, we see explicitly that decaying modes decay time- 
algebraically and not exponentially, in agreement with the lack of spec- 
tral gap. Moreover, the low-frequency/large-time-space behavior is approxi- 
mately hyperbolic, propagating along characteristics associated with the end- 
states u±. 

1.2.2 Evans/Lopatinski determinants and transition to instability 

Further information may be obtained at the spectral level by the comparison 
of Evans and Lopatinski determinants D and A for the viscous and inviscid 
problems (11.2j) and (11.11) . The Evans function D(\), defined as a Wronskian 
of functions spanning the decaying manifolds of solutions of the eigenvalue 
equation 

(L -X)u = 

associated with L at x —>■ +oo and x — > — oo is an analytic function with 
domain containing {9?A > 0}, whose zeroes away from the essential spectrum 
correspond in location and multiplicity with eigenvalues of L with respect to 
any L p . Its behavior is also closely linked with that of the resolvent kernel 
of L, i.e., the Laplace transform with respect to time of the Green function 
G] see jXGJl [Gil [ZHl ESI EB E2] for history and further details. The 
corresponding object for the inviscid linearized problem is the Lopatinski 
determinant A (A) defined in [Kj IMll IM2| IMej . a homogeneous function of 
degree one: in the one-dimensional case, just linear. 

An important relation between these two objects, established in |GZ| IZS| 
MeZj, is the low- frequency expansion 

(1.13) £>(A) = 7 A(A) + o(|A|) 

quantifying the above observation that low-frequency behavior is essentially 
hyperbolic, where 7 is a constant measuring transversality of the profile u 
as a connecting orbit of equilibria u± in the associated traveling-wave ODE, 
nonvanishing for transversal connections. 
A corresponding expansion 

(1.14) ZWA)= 7 A C j(A) + o(|A|) 

holds for the Evans function D rNS associated with detonation wave solutions 
of (11.6)) . where A C j denotes the Lopatinski condition, not for (|1.5p . but 



9 



for the still simpler inviscid-instantaneous-reaction-rate Chapman- J ouguet 
model in which detonations are modeled by piecewise constant solutions 
( 11.31) across which combustion proceeds instantaneously; see |Z1[ LyZl , LyZ2 



IJLWj ICJLW] for further details. That is, the low-frequency behavior is not 
only hyperbolic but "instantaneous", with small-scale details of the profile 
structure lost. Both (I1.13P and (1 1.141) extend to the corresponding mul- 
tidimensional problem of a planar shock moving in M. d (different from the 
finite-cross-sectional case considered here). 

As calculated in |Mlj and |LyZl| , respectively, neither A nor Acj vanishes 
for A 7^ for an ideal-gas equation of stateQ Likewise, 7 does not vanish 
for any choice of parameters for ideal gas dynamics by a well-known result 
of Gilbarg [GI], or, by results of [GSJ for reactive gas dynamics with ideal 
gas equation of state in the ZND-limit \B\ — > 0. From these observations, 
combined with fll.l3l) - fll.l4l) . we may deduce that D and D r ^s have for all 
choices of physical parameters precisely one zero at A = 0, corresponding to 
translation-invariance of the background equations. 

In particular, as physical parameters are varied, starting from a stable 
viscous traveling-wave, transition to instability, signalled by passage from 
the stable complex half-plane {3?A < 0} to the unstable complex half-plane 
{3?A > 0} of one or more eigenvalues of the linearized operator about the 
wave, cannot occur through passage of a real eigenvalue through the origin, 
but rather must occur through the passage of one or more nonzero complex 
conjugate pairs through the imaginary axis: that is, a Poincare-Hopf-type 
configuration. This observation, made in |LyZl , LyZ2| , gives rigorous cor- 
roboration at the spectral level of the numerical observations of Kasimov and 
Stewart [KSJ. What remains is to convert this spectral information into a 
rigorous nonlinear existence result. 



1.2.3 Formulation of the problem 

From physical/numerical observations, one expects for an ideal gas equation 
of state that such transition seldom or never occurs for shocks but frequently 
occurs for detonations. However, from a mathematical point of view, the 
situation for shocks and detonations appears to be entirely parallel. We may 
thus phrase a common mathematical problem: 

1 Indeed, this holds for multidimensions as well; see [Mil IJLW] , 
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(P) Let u £ denote a family of traveling-wave (either shock or detonation) 
solutions indexed by bifurcation parameter e G K 1 , with associated linearized 
operators L(e) transitioning at e = from stability to instability as follows. 

Assuming the "generic" spectral situation that 

(i) the L 2 essential spectrum of each L(e) is contained in {9?A < 0} U {0} ; 

(ii) the translational zero-eigenvalue of each L(e) is simple in the sense that 
the associated Evans function D e vanishes at A = with multiplicity 
one, and 

(iii) a single complex conjugate pair of L 2 eigenvalues of L(e) (zeroes of D e ) 

A±(e) = 7 (e) ± ir(e), 7 (0) = 0, r(0) ^ 0, (dj/de)(Q) > 

crosses the imaginary axis with positive speed, and the rest of the L 2 
point spectrum of L(e) is located in {-RA < 0}, 

show that there occurs a Poincare-Hopf bifurcation from the family of trav- 
eling wave solutions u e to nearby time- oscillatory solutions. 

We examine this problem in a series of different contexts. 
1.3 Model analysis I: the scalar case 

Let us first recall the situation considered in [TZ1] , of a one-parameter family 
of standing- wave solutions u 6 (x) of a smoothly- varying family of equations 

(1.15) u t = F(e, u) := u xx - F(e, u, u x ) 

(possibly shifts F(e,u,u x ) := f(u,u x ) — s(e)u x of a single equation written 
in coordinates x — > x — s(e)t moving with traveling-wave solutions of vary- 
ing speeds s(e)), with linearized operators L(e) := dT jdu\ u= ^ for which 
a spectral gap may be recovered in an appropriate exponentially-weighted 
norm 

WfWk ■= E \\{d/dxyf{x)f LV \\f\\ L > := \\e^ + ^ 1/2 f(x)\\ L2 , V > 0. 

3=0 

We call this the weighted norm condition. 
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This approach, introduced by Sattinger [S] applies to the case (see Sec- 
tion 11.2.11) that signals in the far field are convected under the linearized 
evolution equation Ut = L(e)u inward toward the background profile, hence 
time-exponentially decaying in the weighted norm || • ||#2 penalizing distance 
from the origin. The method encompasses both (11.21) and (jl.6p in the scalar 
case u, z G M 1 , with artificial viscosity B = D = I. It has also interesting 
applications to certain reaction-diffusion systems and the related Poincare- 
Hopf phenomenon of "breathers" [NM, IN} IIIMj . However, it does not apply 
to either shock or detonation waves in the system case u G M n , n > 2. 

In this case, we have the following rather complete result, including sta- 
bility along with bifurcation description. 

Proposition 1.1 ([TZl]). Let u e , ( I1.15P be a family of traveling-waves and 
systems satisfying the weighted norm condition and assumptions (P), with 
F G C 4 . Then, for a > sufficiently small and C > sufficiently large, 
there are C 1 functions e(a), e(0) = 0, and T*(a), T*(0) = 2it/t(0), and a 
C 1 family of solutions 



of (I1.15P with e = e(a), where u a (-,t) is time-periodic with period T*(a) and 
a a is a constant drift, such that 



for all t > 0. Up to fixed translations in x, t, for e sufficiently small, these 
are the only nearby solutions of this form, as measured in H^. Moreover, 
solutions u a are time- exponentially phase- asymptotically orbitally stable with 
respect to if de/da > 0, in the sense that perturbed solutions converge 
time- exponentially to a specific shift in x and t of the original solution, and 
unstable if de/da < 0@ 

Proposition 11.11 was established in |TZlj by center-manifold reduction, 
with the main issue being to accomodate the underlying group invariance of 
translation. The basic idea is to coordinatize u as (v,a), where a param- 
eterizes the group invariance, then work on the quotient space v, reducing 
the problem from relative to standard Poincare-Hopf bifurcation, afterward 

2 In (|1.17p . we have repaired an obvious error in |TZlj . where u° appears in place of u e , 
and in (|1.16p eliminated a redundant time-periodic translation 8 a (t) . 
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recovering the location a (a function of a and t) driven by the solution on 
the quotient space, by quadrature. The integral of the periodic driving term 
yields a periodic part that may be subsumed in the profile and a drift a a t. 
See |TZ1] for further discussion and details. 

Remark 1.2. A consequence of (I1.17P (by Sobolev embedding) is 

(1.18) \u a (x,t) -u £ \<Ce~ vlxl . 

That is, the existence result in weighted norm space includes quite strong 
information on the structure of the wave. On the other hand, the stability 
result is somewhat weakened by the appearance of a spatial weight, being 
restricted to exponentially decaying perturbations. 

1.4 Model analysis II: systems of conservation laws 

The purpose of the present paper is to extend the analysis initiated in |TZ1] 
for scalar models to the more physically realistic system case. For simplicity, 
we restrict to the somewhat simpler case of viscous shock solutions of systems 
with artificial viscosity B = I: however, the method of analysis in principle 
applies also in the general case; see discussion, Section 11.71 

Specifically, consider a one-parameter family of standing viscous shock 
solutions 

(1.19) u(x,t) = u £ (x), lim u £ (z) = u% (constant for fixed e), 

z— >±oo 

of a smoothly- varying family of conservation laws 

(1.20) ut = Tie, u) := u xx - F{e, u) x , u E W 1 , 
with associated linearized operators 

(1.21) He) := -f\ u=a c = -d x A £ (x) + &l 

ou 

denoting 

(1.22) A £ (x) := F u (u £ (x),e), A £ ± := lim A £ (z) = F u (u £ ± ,e). 

We take u £ to be of standard Lax type, meaning that the hyperbolic 
convection matrices A £ + and A £ _ at plus and minus spatial infinity have, 
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respectively, p—1 negative and n — p positive real eigenvalues for 1 < p < n, 
where p is the characteristic family associated with the shock: in other words, 
there are precisely n — 1 outgoing hyperbolic characteristics in the far field. 

This is the only type occurring for gas dynamics with standard (e.g., ideal 
gas) equation of state; for reacting gas dynamics, the corresponding object 
is a strong detonation, which is the only (nondegenerate) type occurring in 
the ZND limit B, D — > [GS] . The special features of Lax-type shocks (resp. 
strong detonations), as compared to more general undercompressive shocks 
(resp. weak detonations) that can occur in other settings, turn out to be 
important for the analysis (see Remark 13.101) . in sharp contrast with the 
generality of |TZlj . 

Note, for the system case n > 2, that there is at least one outgoing 
characteristic, so that the weighted norm methods of the previous section do 
not apply. In particular, we see no way to construct a center manifold for 
this problem, and suspect that one may not exist; the slow (time-algebraic) 
decay rate of outgoing modes evident in our description of the Green function 
in Section 11.2.11 suggests an essential obstacle to such construction. 

On the other hand, the conservative form of system (ll.20p implies that 
relative mass j u(x)dx is conserved for all time for perturbations u = u — u £ , 
u satisfying (jl.20p . Thus, at least at a formal level, there is a convenient 
invariant subspace of (11.201) consisting of perturbations with zero excess 
mass, on which which the zero eigenvalue associated with translational in- 
variance is removed. For, recall that j u'(x)dx = u + — w_ 7^ 0, so that 
the associated zero-eigenfunction is not in the subspace. On the other hand, 
nonzero eigenfunctions always have zero mass [ZHJ, since \<f) = L(e)(f> implies 
A J (p(x)dx = by divergence form of L(e), hence the crossing nonzero imag- 
inary eigenvalues A±(e) in (P) persist. Thus, the spectral scenario (P) trans- 
lates in the zero-mass subspace to a standard Poincare-Hopf scenario, with 
no additional zero-eigenvalue, for which the translational group-invariance 
need not be taken into account. Recall, that this was the main issue in the 
analysis of the scalar case. 

In short, the two problems (scalar vs. system) have essentially comple- 
mentary mathematical difficulties, hence little technical contact. Accordingly, 
the analysis has a quite different flavor in the system case, depending on both 
the full, pointwise Green function bounds of |MaZ3j and the special, conser- 
vative structure of the equations. 

Our result in this case, and the main result of the paper, is as follows. 
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Theorem 1.3. Let u e , (11.201) be a family of traveling-waves and systems sat- 
isfying assumptions (P), with F £ C 2 . Assume further that u £ is of Lax type, 
with o~(A±) real, nonzero, and simple, A± as defined in ( 11.221) . Then, for 
a > sufficiently small and C > sufficiently large, there are C 1 functions 
e(a), e(0) = 0, and T*(a), T*(0) = 2tt/t(0), and a C 1 family of solutions 
u a (x,t) of ( ll.20p with e = e(a), time-periodic of period T*(a), such that 

(1.23) G~ x a < sup(l + \x\)\u a {x,t) -w e(a) (x)| < Ca, forallt>0. 

Up to fixed translations in x, t, for e sufficiently small, these are the only 
nearby solutions as measured in norm \\f\\xi '■= 11(1 + I^D/Wllz/ 00 ^) that are 
time-periodic with period T £ [T ,Ti], for any fixed < T < T\ < +oo. If 
u e + 7^ u e _, they are the only nearby solutions of the more general form ( 11.161) . 

Note that the statement of Theorem 11.31 is considerably weaker than that 
of Proposition II. 1[ asserting no stability information, and only the relatively 
weak structural information of algebraic decay (I1.23P at ±oo, as compared 
to the exponential bound (ll.lSp . As suggested by our formal "zero-mass" 
discussion, the periodic solutions constructed have zero mean drift a, in 
contrast to (I1.16p . However, somewhat surprisingly, they do not appear to 
have zero excess mass. See Remarks 11.81 and 12.101 for further discussion and 
explanation. 

1.5 Model analysis III: flow in an infinite cylinder 

One-dimensional traveling-wave solutions (11.41) of (11.21) with B = I may 
alternatively be viewed as the restriction to one dimension of a planar viscous 
shock solution 

(1.24) u(x,t) =u(x 1 - st) 

of a multidimensional system of viscous conservation laws 

(1.25) u t + Y^ f J (")■<-, = &xU, u £ IR n , x £ M. d , t £ R + 

on the whole space. Likewise, traveling-wave solutions (ll.24p may be viewed 
as planar traveling-wave solutions of (11.251) on an infinite cylinder 

C := {x : (xi, x) £ M, 1 x Q}, x = (x 2 , . . . , Xd) 
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with bounded cross-section G K , under artificial Neumann boundary 
conditions 

du/dx ■ v n = for x G <9fi, 

or, in the case that Q is rectangular, periodic boundary conditions, Q = T d_1 , 
T d ~ l the rectangular torus. 

We take this as a simplified mathematical model for flow in a duct, in 
which we have neglected boundary-layer phenomena along the wall dQ in 
order to isolate the oscillatory phenomena of our main interest. Consider a 
one-parameter family of standing planar viscous shock solutions u £ (xi) of a 
smoothly-varying family of conservation laws 

d 

(1.26) u t = T[e, u) := A x u - ^ F j {e, u) Xj , u G W n 

3=1 

in a fixed cylinder C, with Neumann (resp. periodic) boundary conditions- 
typically, shifts F J (£, w)^. := ^2P{u) X:j — s(e)u Xl of a single equation 
( I1.25P written in coordinates X\ — > x\ — s(e)t moving with traveling- wave so- 
lutions of varying speeds s(e)- with linearized operators L{e) := dT j du\ u= ^ ■ 
As in the previous subsection, we take u £ to be of standard Lax type, con- 
sidered as a shock wave in one dimension. For simplicity, we take Q = T^ 1 . 
Then, we have the following result generalizing Theorem 11.31 

Theorem 1.4. Let u e , (11.261) be a family of traveling-waves and systems 
satisfying assumptions (P), with F G C 2 , Q = T d_1 . Assume further that 
u e is of Lax type, with o~(A±) real, nonzero, and simple, A± as defined in 
(ll.22p . Then, for a > sufficiently small and C > sufficiently large, there 
are C 1 functions e(a), e(0) = 0, and T*(a), T*(0) = 2?r/r(0), and a C 1 
family of solutions u a (xi,t) of (11.201) with e = e(a), time-periodic of period 
T*(a), such that 

C~ l a< sup (1 + \x x \)\u a {x x ,t) -w £(a) (xi)| < Ca, forallt>0. 

Up to fixed translations in x, t, for e sufficiently small, these are the only 
nearby solutions as measured in norm := ||(l + l a; i|)/( a; )||L oo (a;) that are 

time-periodic with period T G [Tq,Ti], for any fixed < T < T\ < +oo. If 
u e + ^ u e _, they are the only nearby solutions of the more general form (11.161) . 

Remark 1.5. Extension to a circular cross-section Q := {\x\ < R} follows 
by Bessel function expansion as in [KSj . The extension to general cross- 
sectional geometries is discussed in Section 16.41 
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1.6 Idea of the proof 



We now briefly discuss the ideas behind the proof of Theorems 11.31 and 11.41 
Under the spectral assumptions (P) (stated in Section 11.2.31) . there exist 
smooth (in e) L(e:)-invariant projections onto the two-dimensional eigenspace 
S £ of L{e) associated with the pair of crossing eigenvalues A±(e) = 7(e)±ir(e) 
and its complement S £ . Projecting onto these subspaces, and rewriting in 
polar coordinates the flow on E e , we may thus express (11.201) in standard 
fashion as 

r = j(e)r + N r (e, r, 9, v), 

(1.27) 9 = T(e)+N e (£,r,e,v), 

v = L(e)v + N v (e, r, 9, v), 

where the "transverse linearized operator" L(e) is the restriction of L(e) to 
£ e , and Nj are higher-order terms coming from the nonlinear part of fl 1 . 2 p : 
N r and N v quadratic order in r, v and their derivatives, and Ng linear order. 

We are precisely interested in the case that L has no spectral gap, i.e., 
9?(x(L) ^ —9 < for any 9 > 0. One may think for example to the finite- 
dimensional case that L(0) has additional pure- imaginary spectra besides 
A±(0) = ±2t(0). Thus, the center manifold may be of higher dimension, 
and one cannot follow the usual course of reducing to a center manifold 
involving only r, 9, e and applying the standard, two-dimensional Poincare- 
Hopf Theorem. Instead, we proceed by a direct analysis, combining the 
two-dimensional Poincare return map construction with Lyapunov-Schmidt 
reduction. 

1.6.1 Return map construction 

Specifically, truncating \v\ < Cr, C ^> 1, in the arguments of N r , Ng, we 
obtain 

\N \ < C 2 r, 

and therefore 9 > r(0)/2 > for e, r sufficiently small. Thus, in seeking 
periodic solutions 

(1.28) (r,MCO = (r,M(0), 

we may eliminate 9, solving for T(e, a, b) as a function of initial data 

(a, b) := {r,v){0) 
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and the bifurcation parameter e, with T(0, 0, 0) = 2ir/ r(0), and seek solutions 
(I1.28P as fixed points 

(1.29) (a,b) = (r,v)(T(e,a,b)) 

of the Poincare return map (r,v)(T(e, ■, ■)). 

Using Duhamel's formula/variation of constants, we may express (11.291) 
in a standard way (see, e.g., [HKj ) as 

= f(e, a, b) = ( e ~f(c) T (e,a,b) _ 1 \ a + N 

(1 30) 

= g(e, a, 6) = ( e ^)r(*.«.*) - Id) 6 + JV 2 , 



where 



(1.31) 



<-T(e,a,b) 

Nf.= / e^ £){T{£ > a > b) - s) N r (s,r,6,v)(s)ds, 



o 



rT{e,a,b) 

N 2 := / e L(£){T(e ' a ' b) - s) ^(e,r,^,t;)(s)ds 



o 



are, formally, quadratic order terms in a, b. Note that (r, 9, v) are functions 
of (e, a, b) as well as s, through the flow of fll.27j) . with 

||M)||<C||(a,&)|| 

for any "reasonable" norm in the sense that (I1.27P are locally well-posed (in 
practice, no restriction). By quadratic dependence of Nj, we have, evidently, 
that (e, a, b) = (e, 0, 0) is a solution of fll.30p for all e. 



1.6.2 Standard reduction 

Continuing in this standard fashion, we should next like to perform a Lyapunov- 
Schmidt reduction, using the Implicit Function Theorem to solve the g equa- 
tion for b in terms of (e, a), i.e., to find a function 

b = B(e,a), £(0,0) = 0, 

satisfying g(e, a, B(e, a)) = 0: equivalently, to reduce to the nullcline of g. 

If this were possible, from |A^ 2 | < C\(a,b)\ 2 in (11.30D (ii). we would find, 
further, that \B(e, a)\ < C3|a| 2 , from which straightforward Duhamel/Gronwall 
estimates on (11.271) (i).(ii) would yield 

\v\ < C 2 r 2 
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for s € [0, T], justifying a posteriori the truncation |u| < Cr, of Ng performed 
in the first step, provided that (e, a, b) are taken sufficiently small: specifi- 
cally, small enough that r remains less than or equal to C/C2 for s G [0, T]. 

Substituting into the / equation, we would then obtain a reduced, scalar 
bifurcation problem 

= /*(£, a) := /(e, a, B(e, a)), f*(e, 0) = 

that would be solvable in the usual way (i.e., by dividing out a and applying 
the Implicit Function Theorem a second time, using 

d £ (a-\)(0,0) = d e (j/r)(0) = <9 e7 (0)/t(0) ^ 

to solve for e as a function of a; see, e.g., [HK] . or Section I2.3.ip . The 
Poincare return map system (11.301) would then be solved, and a solution of 
(ll.30p would generate a periodic solution of (11.271) . 

1.6.3 Regime of validity 

Let us ask ourselves now under what circumstances this basic reduction pro- 
cedure may actually be carried out. Formal differentiation of (I1.30P yields 

d b g(0, 0, 0) = e 27r(Z/r)(0) - Id. 

In the finite-dimensional ODE case, the usual condition of application of the 
(standard) Implicit Function Theorem is thus that (e 27r ( L//T ^°) — Id) be invert- 
ible, i.e., that L(0) have no "resonant" oscillatory modes, i.e., pure imaginary 
eigenvalues that are integer multiples of the crossing modes A±(0) = ±ir(0). 
Note that this includes interesting cases for which L has no spectral gap, 
and thus that do not yield to the standard center-manifold reduction. The 
analogous criterion in the infinite-dimensional (e.g., PDE) setting is that all 
possible resonant modes mr(0), nGZ, lie in the resolvent set of L. 

This is encouraging, and shows that the requirements of Lyapunov-Schmidt 
reduction are much less than those for center-manifold reduction: in partic- 
ular, only spectral separation (from both ±ir(0) and their aliases mr(0)) 
rather than spectral gap is needed. Unfortunately, in the setting (ll.20p . 
(I1.26P of our interest we have not even spectral separation, since essential 
spectra of L(e) accumulates at A = for every e. In this case, or others for 
which separation fails, we must follow a different approach. 
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1.6.4 Refined analysis 

Alternatively, we may rewrite (11.301) (ii) formally as a fixed-point equation 

(1.32) b= (ld-e £ ^ T ( £ '°' 6 )) -1 ^ 2 (£,a,&), 

then apply the Contraction-mapping Principle to carry out an Implicit Func- 
tion construction "by hand" (reminiscent of, but not exactly the standard 
proof of the Implicit Function Theorem). In the standard case that \e Lt \ is ex- 
ponentially decaying (~ spectral gap), we may rewrite (11.321) using Neumann 
expansion as 

oo 

(1.33) b = e^ e)T{£ ' a ' b) N 2 {e, a, b). 

j=0 

The basis for our analysis is the simple observation that, more generally, 
provided the series on the righthand side converges (conditionally) for || (a, b) \\ 
bounded, then (i) (by scaling argument) the righthand side is contractive 
in b for ||(e, a, b)\\ sufficiently small, and (ii) (by standard telescoping sum 
argument, applying (id — e L ^ T ^ £,a,b ^ to (jl.32p ) the resulting solution 

b = B(e,a) 

guaranteed by the Contraction-mapping Principle is in fact a solution of the 
original equation fll.30l) (ii). See Section [2] for further details. 

In the context of (ll.20p . by divergence form of the equation, N 2 = (n 2 ) x , 
with n 2 = 0(\(r,v)\ 2 ). For simplicity of discussion, model n 2 as a bilinear 
form in (r, v), so that (by quadratic scaling) n 2 G L 1 for v G L 2 , from which 
we obtain N 2 = n x , n G L l for b G L 2 . Thus, replacing the righthand side of 
(11.33P by its continuous approximant 

/'OC /"OO 

N 2 + T l e L{£)t N 2 dt = N 2 + T l {e i{£)t d x )ndt, neL 1 , 

JT JT 

we see that convergence (as a map from L 2 to L 2 ) reduces, roughly, to a 
space-time stability estimate 

/ G y (x,t;y)dt)n(y)dy < C\\n\\ L i 

-oo K JT J L 0*0 
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where G is the Green kernel associated with transverse solution operator e Lt . 
Estimate (1 1.341) is quite similar to those arising for the full solution operator 
(again, through Duhamel's principle; see, for example, [Z4, MaZ4j) in the 
study of nonlinear stability of spectrally stable waves. That (11.341) holds for 
all n 6 L 1 , is equivalent to 



POO 

1-35) sup / G y (-,t;y)dt 

y Jt 



< c. 

LP- 



(One implication is a consequence of the triangle inequality; for the other 
take a sequence converging to a Dirac function in L 1 .) This would hold, for 
example, if 



POO 

1.36) sup / \\G y (-,t;y)dt\\ L2 < C. 

y Jt 



Note that the neglected first term N2, corresponding to j — in (ll.33p . 
is C 2 from L 2 to L 2 by the smoothing action of (ll.3ip (ii) ; see Section HI 

In the arguments above, we have used strongly the specific structure 
of the nonlinearity, both_ quadratic dependence and divergence form, since 
the L 2 -operator norm \e Lt \ L2 _ fL 2 does not decay in t. However, this is still 
not enough. For, recall the one- dimensional discussion of Section 11.2. 1} 
modeling G qualitatively as the sum of a convected heat kernel K as in 
f ll.lOP and an error-function term J as in (11.111) . From the standard bounds 
\\K y (;t;y)\\ L 2 = Ct~ 3/4 , \\J y (-,t;y)\\ L 2 < Ct' 112 , hence (the support of K 
and J being essentially disjoint) C° G y (x,t;y)dt is not absolutely conver- 
gent in L 2 (x), and the bound (11.361) does not hold. 

Instead, we must show conditional convergence, using detailed knowledge 
of the propagator G to identify cancellation in G y (x, t; y)dt. For example, 
using K y = a^i^Kt — K yy ), we find that 



;i.37) 



/ K y (x,s;y)ds = a x ( / K t (x, s;y) ds - / K yy (x, s;y) ds) 
Jt Jt Jt 

= a- l {K(x, S]y)\ l T - / K yy (x,t;y) ds) 

Jt 



is the sum of —a^ 1 K(x,T;y) E L 2 {x) with terms a~ 1 K(x,t;y) ~ t~ 1//4 and 
a -1 K yy (x,t;y) ds ~ j/ r s~ 5 ^ 4: ds that are respectively decaying and abso- 
lutely convergent in L 2 (x), hence converges uniformly with respect to y in 
L 2 (x). 
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A similar cancellation argument yields convergence of the error-function 
term, along with uniform boundedness with respect to y. Convergence of 
the error-function term is not uniform with respect to y, a detail which 
necessitates further technicalities: in particular, the weighted estimate 

\v{x)\ < C(l + Ixl)- 1 

leading eventually to (jl.23p . Likewise, there are further issues associated with 
e-regularity of the solution, needed for bifurcation analysis of the reduced 
equation. However, the main idea is contained in calculation (11.371) . 

The mult i- dimensional case goes similarly, with the computation of the 
critical neutral, zero transverse wave number reducing to the one-dimensional 
case, and all others to the case of a spectral gap. 

Remark 1.6. It is readily checked that the above arguments go through in 
the general case that ri2 is quadratic order for \v\l°° < C and not bilinear 
in (v,r), substituting v G L 2 fl L°° (=> 112 G L 1 ) for v G L 2 , working in the 
L 2 n L°° norm in place of L 2 , and substituting for (I1.35P the estimate 



sup 

v 



G y (x, t] y)dt 



< c. 

L 2 C\L°°{x) 



Similar cancellation estimates are important in the study of asymptotic be- 
havior of stable viscous shock waves |L2l ISX} IL3j |R| IHRZj . 

Remark 1.7. As in the result of Proposition ll.ll obtained by center manifold 
reduction, the family of periodic solutions obtained, up to translation, lies 
tangent to the (r, #)-plane corresponding to linearized oscillatory behavior. 

Remark 1.8. As noted earlier, operator S = e LT by divergence form of 
the equations preserves the property of zero mass, J udx = 0, hence each 
finite approximant ^2f =0 to (J — S)^ 1 also preserves zero mass. However, 
f u dx = is not closed with respect to || • \\ L 2, as may be seen by the example 

K(x, S;2 /)|^-K(x,T ;2 /) 

in L 2 as t — ► 00; that is, mass may "escape at infinity" in the limiting process. 
Comparing with (11.371) . we see that the principal term in linijv^oo ^^L SW X 
is exactly of this form, hence the L 2 -limit (J — S)^ 1 does not preserve zero 
mass, even though each finite approximant does. Accordingly, there is no 
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reason that the solution b = b(a, e) = (I — S) -1 ^^, a, e) obtained by reduc- 
tion should have zero mass: in the model case n = u 2 , G = K, it can be seen 
that it does not. Thus, the idea of carrying out a simplified bifurcation anal- 
ysis in the invariant subspace of zero-mass perturbations, though appealing 
at a formal level, in general cannot succeed. See also Remark 12.101 



1.7 Discussion and open problems 



Theorems II .31 and 11.41 together with the spectral observations of |LyZl LyZ2| 



give rigorous validation in a simplified context of the formal and numerical 
observations of [BMR} IKSj . An interesting problem for future work is to 
extend these results to the originally-motivating case of detonation waves of 
the full, reacting compressible Navier-Stokes equations. 

We expect that our one-dimensional analysis will extend in straightfor- 
ward fashion, combining tools developed in |MaZ3t IMaZ4l IZlj and jLRTZj 
to treat, respectively, nonreacting gas dynamics with physical, partial viscos- 
ity and reacting gas dynamics with artificial viscosity. Likewise, we expect 
that we can readily treat flow in a cylinder for the physical equations with 
artificial Neumann or periodic boundary conditions. However, the treat- 
ment of physical, no-slip boundary conditions (presumably associated with 
characteristic viscous boundary layers) involves technical and philosophical 
difficulties beyond the scope of the present analysis, as yet unresolved even 
for nonreacting, incompressible flows: for example, even the construction of 
a background traveling profile becomes problematic in this setting. We point 
out that viscous boundary effects (since also viscosity) are neglected also in 
the ZND setting of jBMRl IKS] . 

A second natural direction for future investigation is the question of sta- 
bility of the periodic waves whose existence we have established here. In 
the absence of a spectral gap, our method of analysis does not directly yield 
stability as in the case of center manifold reduction, but at best partial 
information on the location of point spectrum associated with oscillatory 
modes, with stability presumably corresponding to the standard condition 
de/da > 0. The hope is that we could combine such information with an 
analysis like that carried out for stationary waves in |ZHl IMaZ3j . adapted 
from the autonomous to the time-periodic setting: that is, a generalized Flo- 
quet analysis in the PDE setting and in the absence of a spectral gap. We 
consider this a quite exciting direction for further development of the theory. 

Bifurcation in the absence of a spectral gap has been considered by a 
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number of authors in different settings; in particular, it has been studied 
systematically by loss et al HH H21 HXJ US HM1 ESD ES2] in various contexts 
using an alternative "spatial dynamics" approach. It would be very inter- 
esting to investigate what results could be obtained by this technique in the 
context of viscous shock waves: more generally, to relate it at a technical 
level to the one used here. Considered in this larger context, the interest of 
the present approach is that it gives a simple and explicit connection between 
stability and bifurcation, in the spirit of center-manifold reduction and for- 
mal asymptotics, but adapted to the boundary case of zero spectral gap and 
time-algebraic decay. On the other hand, the price of this "direct" approach 
is that the needed estimates may in practice (as here) be rather delicate to 
obtain. 

We mention also a recent work of Kunze and Schneider [KuS] in which 
they analyze pitchfork bifurcation in the absence of a spectral gap using Sat- 
tinger's weighted- norm method, as described in Section fl~3| but in a situation 
where convection is outward, away from the profile layer. This entails the 
use of "wrong-way" exponentially decaying weights at the linearized level, 
introducing a spatially-exponentially growing multiplier in quadratic-order 
source terms, in combination with separate, compensating estimates at the 
nonlinear level. For similar arguments in the context of nonlinear stability, 
see, e.g., \PW\ IDoj . This interesting approach has been used successfully 
in the shock- wave context to treat the scalar undercompressive case |Doj ; 
however, it does not appear to generalize to the system case. 

Remarks 1.9. 1. Birtea et al in |BPRT] successfully carry out a rather 
complicated bifurcation analysis in the ODE setting without explicit knowl- 
edge of the background solution. Similar techniques might perhaps be useful 
in treating flow in a duct with physcial, nonslip boundary conditions, for 
which description of the background flow is itself problematic. 

2. In carrying out our nonstandard Implicit Function construction for 
shock waves, we faced the problem of non-uniform convergence of series ( 11.331) 
due to lack of spatial localization. We remedied this problem by additional 
weighted-norm estimates. However, the example of another famous non- 
standard Implicit Function construction, namely, Nash-Moser iteration, sug- 
gests the alternative approach of introducing a "localizing" step dual to the 
smoothing step in the Nash-Moser scheme. It would be interesting to see if 
this approach could also be carried out, thus avoiding the need for additional 
analysis associated with pointwise bounds. We note that Nash-Moser iter- 
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ation combining temporal with the usual frequency cutoffs has been carried 
out by Klainerman [Klj in the context of a nonlinear wave equation. 

Plan of the paper. In Section [2J we formalize the reduction procedure 
set out in Section 11.6.41 as an abstract bifurcation framework suitable for 
application to general discrete dynamical systems in the absence of spectral 
gap. In Section [3X we recall the pointwise estimates furnished by the meth- 
ods of |ZH[ IMaZ3j on the Green function G associated with the transverse 
linearized solution operator e Lt , then, in Sections @] and use these to verify 
that the shock bifurcation problem indeed fits the hypotheses of our abstract 
framework. This establishes a Lipschitz version of Main Theorem ll.3| we 
improve this to C 1 by a bootstrap argument in Section 15.41 In Section [6l we 
describe the extension to multidimensions, verifying Theorem 11.41 

Acknowledgement. Thanks to Claude Bardos for pointing out the ref- 
erences [HI |I2], and to Bjorn Sandstede and Arnd Scheel for general discus- 
sions of the spatial dynamics method. Thanks to Walter Craig for pointing 
out the reference [Klj . to Thierry Gallay for pointing out [KuS] . and to Tu- 
dor Ratiu for pointing out [BPRT] . K.Z. thanks the Ecole Polytechnique 
Federale de Lausanne for their hospitality in the course of a two-week visit 
within their 2006 special semester on fluid dynamics, during which much of 
the writing was done. B.T. thanks Indiana University for their hospitality 
during several collaborative visits in which the main analysis was carried out. 

Notes. Since the completion of this work, there have been several further 
developments. In |SS3] , Sandstede and Scheel recover and somewhat sharpen 
our results using spatial dynamics techniques, answering the question posed 
in Section [TT71 of what results may be obtained by these methods, obtaining 
the additional information of exponential localization of solutions and ex- 
change of spectral stability. In [TZ31 ITZ3j . we extend our results to shock 
and detonation waves of systems with physical viscosity, at the same time 
greatly sharpening and simplifiying the basic cancellation estimate. The lat- 
ter yields also exponential localization and appears to shed light on the more 
general question of the technical relation between spatial dynamics methods 
and the "temporal dynamics" method used here; see Remark 1.10, [TZ3J. 
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2 Abstract bifurcation framework 



We begin by formalizing the approach described in Section H.6.4[ providing an 
abstract framework for the analysis of bifurcations in the absence of a spec- 
tral gap for discrete dynamical systems in the general form (11.301) : that is, 
the portion of our analysis occurring after the reduction via Poincare return- 
map construction to a fixed-point problem. The specifics of the return-map 
construction are discussed separately in Section HI We carry out the analysis 
at the level of Lipschitz regularity; this can be improved to C r regularity by 
substituting everywhere C r for Lipschitz regularity in the various assump- 
tions. 



2.1 Generalized Lyapunov— Schmidt reduction 

Given two Banach spaces X and Y, we denote by C(X, Y) the space of linear 
continuous functions from X to Y, and by C(X) the space of linear continuous 
functions from X to X. Let B\, B 2 , X\, X 2 be Banach spaces, 

X 2 ^X 1 ^B 1 , x 2 ^B 2 ^B 1 , 

with norms 

II • ||si < || • |Ui < || • IU 2 , II ■ 1 1 Si < II • \\b 2 < II • IU2, 

and for which the unit ball in X% is closed in B\. 
Consider a system of difference equations 

Aa = f(e,a,b) = (R(e,a,b) - Id) a + iVi(e,a,6), 
Ab = g(s, a, b) = (S(e, a, b) - Id)6 + N 2 (e, a, 6), 

associated with a discrete dynamical system 

a = a + f(e, a, b) = R(e, a, b)a + N\{e, a, b), 
b = b + g(e, a, b) = S(e, a, b)b + N 2 (e, a, b), 



(2.2) 
with 



(/,(/)(e,0,0) = 0, 

where e G M m is a bifurcation parameter and R(e, a, b) and S(e, a, b) are 
primary and transverse "linearized" solution operators for one time-step of 
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R : (s,a,b) G R m x R n x X 1 -> R(e,a,b) G £(R n ), 

5: (e,a,6) G R m x R" x X x ^ S(e,a,b) e C(Bi) f] C(X X ), 

and iVi and iV2 are nonlinear maps 

Ni : (e, a, 6) G R m x R" x X 1 — > Ni (e, a, 6) G R n , 
N 2 : (e, a, b) G R m x R™ x X x -> JV 2 (e, a, 6) G X 2 . 

We assume that i?, iVi and _/V 2 are continuous in (e,a,b), as maps from 
R m x R n x B x to £(R n ), R n and X 2 respectively. 

We assume moreover that Ni and N 2 are of quadratic order, and Lipschitz, 
in (a, b), with 



(2.3) 

\d a Ni(e,a,b)\c(B»jB») + |<kNi(e,M)U(Bi,R») < C(|o| + 



2 

i y ' 



|JVi(e,a,6)| < C(|a| + ||6|U 

k), 

and 

(24) ||AT 2 (e,«,6)|U 2 < C(|a| + H&IUJ 2 

\d a N 2 {e,a,b)\ m «,B2) + \dbN2{e,a,b)\c( Bl ,B2) <C(\a\ + \\b\\ Xl )- 



We also assume that R satisfies the bound 

(2.5) \R(e,a,b) - R(e,0,0)\ < C(\a\ + \\b\\ Xl ), 

and that R is Lipschitz with respect to a and b. 

In Sections 12.1.31 and I2.2[ we assume in addition the e-Lipschitz bounds 



(2.6) IdeN^e, a, &)| £(Rm , R ») < C(\a\ + \\b\\ Xl 



) 2 - 



(2.7) \d e N 2 {e,a,b)\c^,B,) < C(\a\ + \\b\\ Xl ) 2 , 
and 

(2.8) \d £ (R(e,a,b) - R(e,0,0))\ < C(\e\ + \a\ + \\b\\ Xl ). 

We assume that the constants in all the above estimates are uniform in 
e,a,b, for |e| + \a\ + \\b\\ Xl small enough. 
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The domains of definition of R, S, N\ and N 2 , together with the domain 
of validity of the above estimates, can be restricted to a neighborhood of the 
origin in R m xR"xIi. 

We have in mind a general discrete dynamical system whose linearized 
solution operator has a finite-dimensional subspace a that is spectrally sep- 
arated from its complement b at e = 0, in which case form (12.21) may always 
be achieved, with R = R(e) and S = S(e). We allow dependence of R, S on 
(a, b) as well as e in order to admit the slightly more general class of systems 
fll.30p derived by the return map construction of Section 11.6.11 

The case of our particular interest is bifurcation from a simple eigenvalue, 
of equilibria of (12. 2p . or, equivalently, zeroes of (12.11) . which corresponds to 
m = n = 1, R(e, 0, 0) differentiable at e — 0, with 

(2.9) R(0, 0, 0) = 1, d £ R(0, 0, 0) ^ 0. 

Remark 2.1. The reason for introducing the different norms X and B is 
that we shall ultimately carry out a fixed-point iteration using the standard 
extension of the Contraction-mapping Principle to the case of a map T that 
is bounded in norm || • \\x 1 and contractive in the weaker norm || • for 
bounded || • ||x i; with the unit ball in X 1 closed in B\. This gives additional 
flexibility that is useful in applications. 

Remark 2.2. In the return map construction of Section 11.6.11 it is more 
standard [HK1 ITZlj to eliminate t instead by using dt/d6 = T(e) _1 + 0(r) 
to rewrite (I1.27P with ""' denoting d/d6 as 

r' = ( 7 /r)(e)r + iVi(r,M,e), 
v' = (L/r)(e)v + N 2 (r,9,v,e), 



leading to a return map (ll.30p of form (12.21) with R, S depending only on e. 



However, the nonlinearity 

N 2 = ((r + Ng)- 1 - t)Lv + (r + Ne^N, = 0{r)Lv + 0(N V ) 

in (12.101) contains terms d x (rv) that do not take v G L 2 to d x L l , to which the 
analysis of Section 11.6.41 does not apply. It is for this technical reason that 
we allow (a, 6)-dependence of R, S, gaining needed flexibility at the expense 
of slight notational and expositional inconvenience. The quasilinear form of 
S may also be useful in more general circumstances, allowing more general 
iteration schemes in situations of limited regularity. 
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Remark 2.3. In the Poincare-Hopf context (see Section H]), e, a G M. and 
R(e,a,b) = e< {£ ^ e ^ with 7 G C k , T Lipschitz, T(0) ^ 0, 7(0) = 0, 
Y(0) 7^ 0, so that (j2T5) . (jZSD and <^M> hold. 



2.1.1 The equilibrium problem 

We seek to solve the equilibrium problem (/, g) = (0, 0) by Lyapunov- 
Schmidt reduction, that is, to determine a function b = B(e,a) satisfying 

g(e,a,B(e,a)) = 0. 

Substituting into (12. ip . we would then obtain a reduced, finite-dimensional 
equilibrium problem (on nullcline, rather than center manifold) 

= f*(e,a) := f(e,a,B(e,a)), 

presumably amenable to analysis by standard finite-dimensional techniques. 

2.1.2 Assumptions 

We are interested in the (non-standard) case that Id — S(e, a, b) does not 
possess a bounded inverse in £(£>i). Our key assumption is as follows. 

Assumption 2.4. For every (e, a, b) G V, a neighborhood of the origin in 
M.™ x M™ x Xi, the linear operator Id — S(e,a,b) G C{Xi) possesses a right 
inverse 

(ld-S(e,a,b))- 1 :X 2 ^B 1 , 
that is bounded in £(B 2 ,Bi) and C(X 2 ,Xi) norms. 

Note that the range of N 2 is contained in X 2 , the subspace of B 2 where 
the right inverse is assumed to be defined. 

Proposition 2.5. For Assumption 12.41 to hold, it is sufficient that, for 
all (e, a, b) G V, the sequence Yl!j=o ^{ £ -> a > be conditionally convergent 
in C(X 2 , Bi) norm, and the limit operator YlJLo S( £ i a i b>Y be bounded in 
C{B 2 , £>i) and C(X 2 , Xi) norms, in which case (Id— S(e, a, b))^ 1 = YlJLo S( £ i a ' 

Proof. Given x G X 2 , let 

N 

(2.11) (I- S^a,^)- 1 ^ := lim V S j {e,a,b)x, 

3=0 
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where the limit is taken in B\. By a standard telescoping sums argument, 
for all N, 

N 

(2.12) (Id - S(e, a,b))J2 S j (e, a, b)x = x - S N+1 {e, a, b)x. 

3=0 

Because S(e,a, b) G £(£>i), 

N N 

(Id — S{e, a, b)) lim N S^{e 1 a ) b)x= lim (Id — S(e, a, b)) y S^{e,a,b)x. 

From the above equality and f!2.12p . 

(Id — S(e, a, &)) (Id — S(e, a, b))~ l x = lim x — S N+1 (e, a, b)x = x, 

as (conditional) convergence of YlJLo ^(e, a, b)x implies that S^(e, a, b)x con- 
verges to 0. Thus, the operator defined by (12.111) is a right inverse of Id — S. 
It is bounded in £(£>2,£>i) and C(X 2 ,Xi) norms by assumption. □ 

We will also need the following regularity assumption. 

Assumption 2.6. For e,a,b in V, the right inverse (Id — S) given by 
Assumption ®. 41 is Lipschitz continuous with respect to b as a function from B\ 
to £(82, Bi) and with respect to (s, a) as a function from M. m+n to £(X 2 , B\), 
in the sense that 

sup ||<9 6 (Id - S') _1 ||£( Bli £( Bs)Bl )) < 00, 

(2.13) 

SUp \\d £ia (ld-S) ||£(Km xR n )iC (x 2) B 1 )) < 00 > 
e,a,b£V 

satisfying also 

sup (I (Id - S)~ 1 (e,a,b)\\ c ^x 2 Jc 1 ) < 00, 

(2.14) £ ' a ' beV 

sup II (Id — S 1 ) (e,a,b)\\c{B2,Bi) < 00. 

£,a,6GV 

Remark 2.7. Convergence in Proposition 12.51 is implied by exponential de- 
cay of 1 1 S n 1 1 c(B2,Bi) as n — > 00, which yields absolute convergence by the ratio 
test. This in turn is implied by stability plus spectral gap of the transverse 
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operator S, i.e., ^(S 1 )! < 1, the standard case. In our case of interest, there 
is no spectral gap and ll^Hr^^), or equivalently ||5 ,n -/V r 2 ||£(B 1 ), is not expo- 
nentially but algebraically decaying, while ||<S' n ||£(,B 1 ) is not decaying at all. 
Thus, convergence will in general be conditional if it occurs, and also depends 
strongly on the specific structure of N 2 . In the application to viscous shock 
waves, moreover, Lipschitz continuity, (I2.13p . holds for the limit, but not the 
approximating finite sums; see Section 15.11 

2.1.3 Reduction 

Lemma 2.8. Under Assumption \2A\ the equation 



is equivalent to 

(2.16) b - (Id - S(e, a, b))- 1 N 2 (e, a, b) £ Ker(Id - S(s, a, b)) n X x . 

Proof. Applying to the left of (12.41) the right inverse (Id — S(e, a, b))~ l given 
by Assumption 12.41 we obtain 

b := (Id - S(e, a, 6))" 1 (Id - S(e, a, b))b = (Id - S(e, a, b))~ l N 2 {e, a, b). 

Observing that b — b belongs to Ker(Id — S(e, a, b)) fl A 1; we obtain (I2.16p . 
Conversely, (12.161) implies g(e, a, b) = 0, by definition of (Id— S(e, a, 6)) _1 . □ 

Consider the map T : R m x R n x X 1 x X 1 ^ X l} defined by, 
T (e, a, 6, 0) = + (Id - S(e, a, 6)) -1 JV 2 (e, a, b). 

Proposition 2.9. Suppose that N 2 satisfies (12.41) and (12. 7p . and that As- 
sumptions^^ and 12.61 are satisfied. Then for \e\ + |a| + sufficiently 
small, T has a unique fixed point b = B(e,a,u) £ X 1} 0,0) = 0, such 
that B is Lipschitz in (e,a,w) with respect to norm \\ ■ with 



(2.15) 



g(e,a,b) = 0, 



(e, a, b) £ R m x R n x X 



(2.17) 



\B(e,a,u)\\ Bl < C{\\u\\ Bl + |a| 2 ) : 
S(e,o,a;)|| Xl <C7(||a;|| Xl + |a| 2 ) 



and, 



(2.18) 



d £ B(e,a,u)\\ c{R ^ Bl) < C(\\u\\ 2 Xi + |a| 2 ). 
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Proof. By (104]) (i) together with (IXD(i). 



11^.0.6.0)11^ ^CdlwII^ + HftH^ + H 2 ) 

for |e| + |a| + |0|x a sufficiently small. Thus, for e, a, fixed and \e\ + \a\ + \0\x 1 
small enough, T maps a small ball 5(0, 2C(|0||x 1 + H)) of Xl centered at 
the origin into itself. 

To prove that T is a contraction, we evaluate — T(b')\\ Bl , for small 

b, b' G X\. This term can be bounded by the sum of 

(2.19) ||(Id-S)- 1 (6)-(Id-5)- 1 (6')||£ (S2 , ei )ll^(6')IU 2 , 
and 

(2.20) \\{ld-S)-\b')\\ c{ B^\N 2 {b)-N 2 {b')\\ B2 . 

By Lipschitz continuity (Assumption I2.4p and (I2.4p (i). for (e,a,b) G V and 
{e,a,b') G V, the term in (I2.19P is bounded by C(\a\ + H&'IUJ 2 !^ - b'\\ Bl , 
where C does not depend on s,a,b,b'. By (I2.14p (ii) and (I2.4p (ii). for \e\ + 
\ a \ + \\b\\xi + H&'Hxi small enough, the term in ( 12.20[) is bounded by C(\a\ + 
\\b\\xi + ||6'||xi)||& — b'Wsx, where C does not depend on e, a, b, b'. 

Combining, we find that T is contractive with respect to b in the B\ 
norm, for (e,a,b) in a neighborhood of the origin in R m x IR n x X\, possi- 
bly smaller than V, and for ||0||xi small enough. By the extension of the 
Contraction-mapping Principle described in Remark 12.11 we may thus con- 
clude the existence of a unique fixed point b = B(e,a,0) G X\ of T in a 
small ball of Xi centered at the origin; specifically, considering the ball in X\ 
as a T-invariant closed set in B\, we may apply the standard Contraction- 
mapping Principle with respect to || • H^. (Recall our assumption that the 
ball in X\ be closed in B\.) The uniqueness of that fixed point, together with 
the identity T(e, 0, 0, 0) = 0, imply that B(e, 0, 0) = 0, for \e\ small. 

By d23D(i)-(ii), we have 

\\B\\ Bl < ||a;|| 01 +C(|a| 2 +|| J B|| 2 1 ), 
\\B\\ Xl < \\to\\ Xl +C(\a\ 2 + \\B\\ 2 Xi ), 

from which, recalling by our fixed point construction that B is confined to a 
small ball in Xi, hence has small B\ and X\ norms, we obtain f)2.17p (i)-(ii). 
Likewise, smallness of B together with (12.141) (ii) gives 

(2.21) sup ||(Id-5 , )- 1 (e,o,S(e,a))|U( ft , Bl ) < oo. 

e,a,beV 
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To establish Lipschitz regularity in e, we evaluate 

II^SII^a, < \\(d e T)(e,B(e))\\ m n Bl) 

+ \\(dbT)(e 1 B(e))\\ c{ B 1 )\\dsB\\ m n A) . 

The second term in the right-hand side of (I2.22p is controlled by contractivity 
of T with respect to b, hence bounded by 0\\d e B\\c( R n tBl ^ for < 9 < 1. The 
first term in the right-hand side of (12.221) is controlled by the sum of 

(2.23) ||(^(W-5)- l )(e,S(e))|| £(B » 1 £ ( xi W )ll^( e » B ( e ))IU« 
and 

(2.24) IICW-^-^e.SCe^ll^^jIK^Ce.SCe))!!^^,). 

By Lipschitz regularity (with respect to e) of (Id — S)^ 1 in C(X 2 , B\) norm, 
(IXTTll(ii). and bound (j23D(i) on JV 2 , the term (|2~23I) is bounded by C(|a| 2 + 
H0IIIJ. Finally, by (12~2H together with the Lipschitz bound (T277D on jV 2 , 
term (I2.24[) is bounded by C(|a| 2 + I^HxJ as well. Combining these bounds 
proves Lipschitz regularity of B with respect to e, with Lipschitz constant 

2C(|a| 2 + \\0\\ 2 Xl )/(l -B). 

Summing all contributions gives (12.181) . Using the facts that (Id — S)~ l and 
N 2 are Lipschitz with respect to a, b, and T is contractive with respect to b, 
we can prove similarly that B is Lipschitz with respect to a and as well. □ 

Noting that G Ker(Id — S(e, a, b)) for any (e, a, 6), we may choose u = 0. 
Then 

g(e,a,B(e,a, 0)) = 0, 

for e, a small. In general, this is the only reasonable choice, since the rest of 
the kernel may depend on the solution b. It follows from Lemma 12.81 that if 
0(e, a) G fl;,Ker(Id — S(e, a, b)), for all e, a, and such that is small enough 
for e, a small, then g(e, a, B(e, a, 0(e, a))) = 0, for all e, a small. 

Remark 2.10. If Id — S has a non-trivial kernel, as holds always in the 
traveling-wave context by translational group invariance, then, in order for 
Id — S to have a right inverse, Range(A*2) must lie in a complementary sub- 
space to Ker(Id — S), which is preserved by S. If also Assumption 12.41 holds, 
then by the choices (I — S)^ 1 = S^, uj = we automatically project out 
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the kernel and work on the (nonlinearly) invariant complementary subspace, 
at least in the case (as for finite-dimensional ODE, or, more generally, when 
zero is an isolated eigenvalue) that the complementary subspace is closed 
with respect to || ■ H^. Note, however, that we do not a priori (i.e., by force) 
restrict to invariant subspaces, but find this in a natural way through the 
analysis. And, in the case that the complementary subspace is not closed 
with respect to || • || J g 1 , we may still obtain a result, even though this result 
may lie outside the complementary subspace. 

For example, in the shock wave context, Range (./V2) is zero-mass, pre- 
served by S, hence (recall, zero-eigenfunction vl has nonzero mass u+ — uJ) 
the set of zero-mass perturbations is a complementary S-invariant subspace 
to Ker(7 — S). However, it is not closed with respect to B\ = || ■ ||x,2, and our 
reduction procedure (apparently) takes us out of this subspace; see Remark 

EH 

Remark 2.11. In the case that S is not stable, it may be split into sta- 
ble/unstable parts and Lyapunov-Schmidt reduction carried out by a com- 
bined forward/backward scheme. In the parabolic PDE context, there can 
exist only finitely many unstable modes, hence the backward unstable flow 
is well-posed and this algorithm may indeed be carried out. 

2.2 Bifurcation 

Extending the discussion at the end of the previous section, consider a Lip- 
schitz (with respect to Xx-norm) curve 0, defined on a neighborhood of the 
origin in R™ x R™, such that 0(e, 0) = and 

(2.25) (e, a) i— » 0(e, a) G p| Ker(Id - ^(e, a, b)) n X u 

\\b\\ Xl <b 

for some b$ > 0. The trivial choice u = shows that at least one such curve 
always exists. 

Proposition 2.12. Under the assumptions of Proposition [279], suppose in 
addition that N± satisfies (12.31) and (12. 6p . that R satisfies (12.51) . (12. 9p and 
( 12.81) . and a,e G R. Then, for each ui as in f |2.25[) . there exists a Lipschitz 
map a i— > 80(a), defined on a ball of radius r centered at the origin in M a , 
such that £0(0) = 0, and 

(/, g)(£0( a )i a , B(e (a), a, 0(e (a), a))) = 0, 
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where B is the Lipschitz map given by Proposition ^ $\ Moreover, for \a\ <r, 
e = 60(a) is the unique solution of (/, g)(-, a, B(-, a, o(-, a)) = (0, 0). Here, the 
size of r depends in particular on the size of the Lipschitz constant of uj. 

Proof. There exists k > 0, such that, for \e\ + \a\ small, ||0||xi < n\a\. 
Consider the map (e, a) 1— ► B := B(e, a, 0(e, a)), where B is given by Propo- 
sition [53B The map B is defined for \e\ + \a\ small, and g(e,a,B) = 0. 
Besides, from (12.171) . the bound on 0, and (I2.3p (i). there holds the bound 

|JVi(e,a,S)| < C{K)\a\ 2 , and 

f(e, a, B) = (R(e, 0, 0) - l)o + (R(e, a, B) - R(e, 0, 0))a + N^e, a, B), 

(recall that #(0,0,0) = 1). Let 

Nx := R(e, a, B) - R(e, 0, 0) + N^e, a, B)/a, N^a = 0) := 0. 

From (I2.5p and the above bound on Ni(B), we see that |iVi| < C\a\, for 
e,a small. The bounds (Q, ([21]) and (T2TT7|) - ff2^8|) imply that N x (B)/a is 
Lipschitz in e, with Lipschitz semi-norm controlled by \a\. The bounds (I2.8p . 
(I2.17P and (I2.18P imply that R(e, a, B) — R(e, 0, 0) is also Lipschitz in e, with 
Lipschitz semi-norm controlled by |e| + |o|. Thus, Ni is Lipschitz in e, with 
Lipschitz semi-norm controlled by \e\ + \a\. Similarly, the bounds (I2.3p (ii). 
(I2.17p - (l2.18p and (12.51) imply that N\ is Lipschitz in a as well. Let now 

/ := r l f{B)/a = e + j- l (R(s, 0, 0) - 1 - js + ^ 

where 7 := <9 e i?(0,0,0) 7^ by (12. 9p . Because R(e,0,0) is assumed to be 
differentiable at the origin, 

(2.26) i2(e,0,0)-l-7£ = o(|e|). 

Thus, for \a\ small enough, the map N := 7~ 1 (-R(e, 0, 0) — 1 — 75 + Ni) maps 
a small ball in to itself. The above description of Ni and (I2.26P imply 
that N is contractive in e, for e, a small enough. The Contraction-mapping 
argument used in Proposition 12.91 then implies the existence and uniqueness 
of a fixed point of —N, in a small ball centered at the origin in R™, for \a\ 
small. This fixed point a 1— > 60(0) is the desired solution to / = 0. Indeed, it 
is Lipschitz in a, because N is, and £0(0) = by uniqueness. □ 

Proposition 12.121 concerns the general case. Let us now consider a more 
specific situation for which Ker(Id — S(e, a, b)) has a uniform structure. 
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Assumption 2.13. The Kernel Ker (Id — S(e, a, b)) is independent ofb, with 
a locally Lipschitz basis V\(e, a), . . . , vg(e, a) of fixed dimension i. 

Under Assumption 12.131 the maps 



are uniformly Lipschitz in (e, a) for a in B(0,R) C M. e and |e| + \a\ small, 
with combined range covering the cone ||u;||xi < C\a\, any C > 0, for R > 
sufficiently large. 

Corollary 2.14. Under the assumptions of Proposition \27[2i suppose further 
that Assumption 12.131 is satisfied. Then, for each \a\ small enough and each 
a G B(0,R) C R , m a neighborhood of the origin in R e x X 1; tnere e:ras£ a 
unique (e, 6) swc/j t/iat 



27ms equilibrium is given by e = e^ia), b = B(e0<r,a,0 a ), with the notations 
of Propositions 12.91 and 12.121 Moreover, all solutions of (f, g) = 
\\b\\x 1 < C|a| and (e, a, 6) sufficiently small are of this type, apart from the 
trivial solution (a, 6) = (0,0). 

Proof. Only the final assertion requires discussion. Evidently, solutions with 
IMIxi < C\a\, u := b— (Id — S(e, a, 6)) _1 A^ 2 (e, a, 6), are of the described type. 
Thus, we have only to observe that ||u;||xi = \\b\\x 1 + 0(\\b\\x 1 + M) 2 , by the 
bounds on N 2 ( f)2.4p (i)). so that ||&||xi < C\a\ implies ||u;||xi < Ci |a| . □ 

That is, there occurs a "restricted" ^-fold Lipschitz bifurcation of equilib- 
ria at e = from the trivial solution (a, 6) = (0, 0), on a cone of expanding 
radius proportional to \a\ around the base solution corresponding to oj = 0. 

2.2.1 Group invariance 

The main situation that we have in mind is the case that Ker(Id — S(e, a, b)) 
corresponds to group-invariance(s) of the underlying equations, so that, in 
particular, the kernel depends only on e. 

Let £ e be the set of equilibria (a, 0) G M™ fl X 1 associated with the value 
e of the bifurcation parameter. Given (a, 6) G £ e , let \\(a, b)\\ := |a| + ||&||xi- 




(f,g)(e,a,b) = 0, 
b - (Id - S(e, a, b))' 1 N 2 (e, a, 6) = 0°(e, a). 
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Assumption 2.15. There exists an £-dimensional Lie Group G, with Lie 
algebra Q, and Eq > 0, such that, for \e\ < Sq, there exists a family of group 
actions 

$ £ : Gxf^ S £ , 

one-to-one in g G G, satisfying the following properties, where the identity 
element in G is denoted by &q: 

(i) The map $ e (-,0,0) is differentiable at ea, uniformly in e, in the sense 
that, for all 7 G C 1 : Q —>■ G, where Q is a neighborhood of the origin 
in R e , 



1 

sup 

|e|<eo 



$ e ( 7 (x), 0, 0) - d^ £ (e G , 0, 0) ■ (f(0) ■ x) 



0, as x —>■ 0. 



(ii) The derivative di$ e (ec,0, 0) is an isomorphism 

G {0r S } x (Ker(Id-5( £ ,a,6))nX 1 ), 

and the norms of ((9i$ £ (ec 0, O))" 1 are bounded uniformly in e, for 

(iii) For g in a neighborhood of ea, the map & 6 {g, •, •) is Lipschitz in (a, b), 
for small ||(a, 6)||, with Lipschitz semi-norm bounded uniformly in e, 
for \e\ < e . 

(iv) For || (a, 6)|| sufficiently small, the map $ £ (-, a, b) is continuous in g for 
g in a neighborhood of e G - 

Proposition 2.16. Suppose that N 2 satisfies (12.41) (i). that Assumptions 12.41 
and 12.151 are satisfied, and moreover that the right inverse satisfies bound 
(I2.14p (i). Then, for any C > 0, there exists p > 0, such that, for each 
equilibrium (a, b) G £ e such that ||(a, fe)|| < p and \e\ < £q, there exists g G G, 
such that (a,b) := & £ (g)(a,b) satisfies \\b\\x 1 < C || a. 1 1 . 

Proof. We first observe, given (a ,b Q ) G £ £ with ||(a ,6o)|| sufficiently small, 
and |e| < Eq, that there exists go such that ||$ e (^o> a, b)\\ minimizes ||$ £ ((7, a, 6)|| 
for g G G lying in a fixed neighborhood jV* of e G - For, by (i)-(ii), ||$ e ((7, 0, 0)|| 
is bounded from below by some cq > for g outside a smaller neighborhood 
M' of e G , whereas, by (iii), \\$ £ (g,a,b) - $ £ (c/,0,0)|| < C||(a,6)||, without 
loss of generality C > 1. Thus, for ||(a,6)|| < c /2C, ||$ £ (^, a, 6)|| > c /2 > 
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|| (a, for g outside A/ 7 , and so, by continuity, (iv), there exists a minimal 

go e M'. 

Let (do, bo) G £ e , and let (a, 6) be such a minimizing equibrium, with 
(a, 6) = $ e ((?o)(ao, M for some go G G. Assume by way of contradiction that 
\ a \ = Ci||6||xi f° r ci > sufficiently small. From g(e,a,b) = 0, by Lemma 
12.81 we have 

b - (Id - 5(e, a, 6))-^ 2 (e, a, 5) = G Ker(Id - o,6))nX 1) 
where, by (£3D(i) and (|2Hjl (i). 

(2.27) ||(Id - ^(e, a, fo))- 1 ^^, a, 6)|U X < C(\a\ + \\b\\ Xl ) 2 , 

uniformly in e. From (I2.27p . we deduce that, if ||o||xi is small enough, 

(a,6) - (0,0) = o(||6|U x ). 

By (ii), associated with is 5 s G Q, such that di$ e (ea, 0, 0) • 5 £ = (0, cj). By 
(i)-(ii), there exists g £ G G, such that 

$ £ ((? e ,0,0)- (0,u) = o(\\u\\ Xl ), 

hence 

$ £ (( 7£ ,0,0)-(a,6)=o(||fe||x 1 ), 

uniformly in e. But, then, Lipschitz continuity, (iii), together with the group 
property, gives 

&{g;\ a, b) = &{g;\ &{g £ , 0, 0) + (||6|UJ) 
= <S>%g;\^(g e ,0,0)) + o(\\b\\x 1 ) 
= o(\\b\\x 1 ), 

uniformly in e, thus contradicting minimality of (a, 6). □ 

Corollary 2.17. Under the assumptions of Proposition I27T21 suppose fur- 
ther that Assumption ®. 151 is satisfied, and that <9i$ e (eG, 0, 0) is Lipschitz in 
e. Then, for (e,a,b) sufficiently small, the solution (eo(a),a,B(eo(a),a,0)) 
determined by = in Proposition 12.121 is the unique solution of (/, g) = 
up to group invariance, apart from the trivial solution (a,b) = (0,0). 
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Proof. Assumption 12.151 implies in particular that Ker(Id — S(e,a,b)) is in- 
dependent of b. If v%, . . . , is a basis of Q, then 



d^ £ (e G , 0, 0) • Vl , . . . , <9!$ £ (e G , 0, 0) ■ v e , 



is a Lipschitz basis of Ker(Id — S(e,a,b)). Thus, Assumption 12.131 is sat- 
isfied, and Corollary 12.141 implies that the solution Eo,b , where b (x) = 
B(eo(a),a,0), and eo is given by = in Proposition 12.121 is unique up 
to group invariance, within the cone ||&||xt < C|a| , for any C > 0. For, by 
equivalence of dimension (Invariance of Domain Theorem), for fixed a, the 
manifold of equilibria constructed in Corollary 12. 14[ indexed by o G M , must 
agree locally with the continuous manifold {Q £ °( a \g, a, B(eo(a), 0))}, indexed 
by g G Mr, since uniqueness implies that the second is contained in the first 
within the cone in question. Finally, Proposition 12.161 states that each equi- 
librium has a representative in the cone \\b\\x 1 < C\a\, giving uniqueness up 
to group invariance with no restriction. □ 

That is, there occurs a complete Mold Lipschitz bifurcation of equilibria 
at e = from the trivial solution (a, b) = (0, 0) in a neighborhood of (e, a, b) = 



Remark 2.18. Note that we carry out a complete relative bifurcation anal- 
ysis without factoring out the underlying group invariance as in [TZlJ. 



2.2.2 Manifold of equilibria 

The ideas of the previous subsection extend readily to the more general case 
that {0} x Ker(Id — S) is the tangent space at zero of an ^-dimensional 
manifold of equilibria, not necessarily generated by group invariance. 
Let || (a, := \a\ + ||6||xi as before. 

Assumption 2.19. There exists an I- dimensional family indexed by 5 G 
of locally invertible Lipschitz coordinate changes $ £,<s : x X\ — >• R" x X\ 
taking ( 12. ip to a system of the same form 



(0,0,0). 



(2.28) 



Aa = (R(e, 6, a, b) — Id) a + Ni(e, 5, a, b), 
Ab = (S(e, 5, a, b) - Id)6 + N 2 (e, 5, a, 6), 



^(8,6,0,0) 



0, such that: 



39 



(i) The maps ^ £ (-) = (* a ,^ 6 )(-) := ($ e '-) _1 (0> 0) and ^ e (-) := ($ E >0(O,O) 
are each differentiable at S = 0, uniformly in e in the sense that 

sup -^-\\^ £ {x) -d s y £ (0) -x)\\ — > 0, as x -»• 0. 

\e\<e \ X \ 

(ii) Ker(Id - S(e, a, b)) = Span-fc^^O)}, while d^(0) is full rank. 

(iii) For 5 in a neighborhood ofO, the map <3> £,<5 (-, •) is Lipschitz in (a, b), for 
small || (a, 

(iv) For \e\, ||(a, 6|| sufficiently small, the map $ £, '(a, b) is continuous in 5 
for \S\ sufficiently small. 

Remark 2.20. Under Assumption l2TW i)-(ii). {0} x Ker(Id - S) is the 
tangent space at 5 = of the manifold of equilibria {$?(8) : 5 G M e }. 

Proposition 2.21. For system (12.281) . suppose that N2 satisfies (I2.4p (i). 
and that Assumptions 12.41 and 12.191 are satisfied, uniformly in S. Then, for 
\e\ < e and any C > 0, for each sufficiently small equilibrium (a, b) of (12.11) . 
there exists S G -6(0, r) such that (a,b) := $ e,<5 (a, b) satisfies \\b\\x! < C||cj||- 

Proof. We first observe, given ||(a ,&o)|| sufficiently small, that there exists 
So such that ||$ e,l5o (a, 6)|| minimizes ||$ e,l5 (a, b)\\ for 5 lying in a neighborhood 
Af about zero. For, by (i)-(ii), H^'^O, 0)|| is bounded from below by some 
Co > for 5 outside a smaller neighborhood Af' of zero, whereas, by (iii), 
$ e,5 (a,&)-$ e '*(0,0) < C||(a,6)||, without loss of generality C > 1. Thus, for 
|| (a, 6)|| < c /2C, \\$ £ > s (a,b)\\ > co/2 > ||(a,6)|| for 5 outside Af 1 , and so, by 
continuity, (iv), there exists a minimal 5q G Af' . 

Let (a , b ) be a sufficiently small equilibrium, and let (a, b) be such a 
minimizing equibrium, with (a, b) = $ e,<s (ao, b ). Assume by way of contra- 
diction that |a| < co||&||xt for Co > sufficiently small, or a = o{b). From 
g(e, a, b) = 0, by Lemma [2.81 we have 

b - (Id - 3(e, a, 6))-^ a (e, a, 6) = u G Ker(Id - S(e, a,b))nXx, 

where, by fl2T4l)(i). 

|| (Id - 5(e, a, ft))" 1 ^^, 0,6)11^ < C(|a| + ||6|| Xl ) 2 = 0(6), 
hence (a, 6) = (0, 0) + 0(0). 
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By (i)-(ii), associated with is 5 G B(0, r) such that ds^f(0) ■ 5 = (0,u>), 
hence 

= (0,0) + o(0) = (a, 6) + 0(0). 
Applying $ e ' 5 to both sides, and recalling the definition of \I/ in (ii), we obtain 

(0, 0) = $ e ^(<5) = $ £ < 5 ((a, 6) + 0(0)) = $ £ - 5 (a, 6) + 0(0), 

or 

$^M) = O (0), 

contradicting the assumed minimality of (a, b). □ 

Corollary 2.22. Under the assumptions of Proposition \2A2\ suppose further 
that Assumption 12.191 is satisfied. Then, for (e, a, 6) sufficiently small, the 
l-dimensional family of solutions (e (a, S), a, B(e (a, 5), a, 0)) determined by 
= in Proposition 12.121 applied to systems (12.281) are i/ie unique solutions 
°f if yd) — 0, apart from the trivial solutions (a,b) G 

Proof. Identical with that of Corollary 12.161 □ 

Remark 2.23. In the context of Hopf bifurcation, existence of such a family 
qe,s amouri t s to the assumption that Ker(Id — S) = Ker(L) corresponds to 
(i.e., is the tangent manifold of) a Lipschitz manifold of equilibria of the 
original differential system, with $ £ ' (5 affine maps consisting of translation 
in u := u + o0 + b by points u e ' 5 on the manifold, composed with linear 
coordinate changes to and from eigenvariables (a,b). 

2.3 Brouwer-based Implicit Function Theorem 

We show in this section that under weaker assumptions, we can still obtain a 
bifurcation result. For Assumption 12.61 we substitute the following, weaker 
assumption. 

Assumption 2.24. The right inverse (Id — S)^ 1 given by Assumption 12.41 
is continuous with respect to e,a,b G V, in £(X2,Bi) norm, and satisfies 
estimates ( 12. 14ft from Assumption 12.61 It is moreover Lipschitz with respect 
to b, as a function from B\ to C(B%, Bi), in the sense that 

(2.29) sup ||9 6 (Id - Sy'WciBudB^)) < 00. 

e,a,beV 
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Lemma 2.25. If S is continuous in (e, a) for (e, a,h) G V with respect 
to the C(X 2 ,Bi) norm and uniformly bounded in C(X 2 ), and the sequence 
^2j =0 S(e, a, b)^ converges uniformly with respect to e, a,b E V in C(X 2 ,B\) 
norm, then (Id — S)^ 1 is continuous with respect to e,a,b E V, in C(X 2 , Bi) 
norm. 

Proof. Since we have assumed Lipschitz continuity in b in the stronger norm 
£(£>i, C(B 2 , £>i), it is sufficient to show continuity in (e, a) for (e, a, b) E V 
with respect to C(X 2 ,B\), which reduces, by continuity of uniform limits of 
continuous functions, to establishing continuity of . Suppressing b, denote 
S' = S{e', a'), S = S(e', a'). Then, 

IKS'?' - si£ ( x 2lBl) < wdsy- 1 - s^s'u^ 

+ \\S 3 ~ (S' - S)\\jC(x 2 ,Bi) 

< \\(sy~ - s , - 7 ~ 1 ||£(x 2 ,Bi)ll<s' / IU(x 2 ) 

from which we obtain the result by induction on j. □ 

The main difference between Assumption 12 . 61 and Assumption 12. 241 is that 
in the latter, we do not assume Lipschitz regularity with respect to e, a, an 
advantage since (Id — S)^ 1 typically depends weakly or not at all on b (see 
Remark 12.21 and Section loTTl) . In place of the Contraction-mapping principle 
of Proposition 12.121 which used Lipschitz regularity in e, a, we will use the 
following Implicit Function Theorem. 

Remark 2.26. In the following fixed point argument (Proposition 12.281) . it 
would be sufficient to assume that for fixed x G X 2 , (Id — S)~ 1 x is continu- 
ous with respect to e, a, b, in B\ norm, which would be implied by uniform 
convergence of ^2f =0 S(e, a, b^x in B\ norm for each fixed x G X 2 . 

Lemma 2.27 (Brouwer-based IFT). Let F : lR m x R n — > R m be continuous 
in a neighborhood of the origin, F(0, 0) = ; with F(S, 0) differentiable at 
5 = 0, and d$F(0, 0) invertible. Then, for \a\ sufficiently small, there exists a 
solution 5 = 5(a) of F '(5(a), a) = 0, 5(0) = 0, with 5(-) continuous at a = 0. 

Proof. Expand 

d 5 F(0, 0)~ 1 F(5, a) =5 + n x (5) + n 2 (5, a), 
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where 



m(«J) := d s F{Q, 0)" 1 (F(6, 0) - F(0, 0) - 8d s F(0, 0)) . 



and 



n 2 {5, a) := d 5 F(0, O)" 1 (F(5, a) - F(5, 0)) . 



Both m and 712 are continuous in a neighborhood of the origin. By differen- 
tiability, rii = o(|5|), while by uniform continuity in a compact neighborhood 
of the origin, n 2 — ■» as |a| — > 0. Thus, for |a| sufficiently small, — ni — 7i2 
takes a small ball centered at the origin in M. m to itself, of radius 



whence there exists a solution 5 = (—ni+n 2 )(5, a) by the Brouwer fixed-point 
Theorem. By construction, 5(a) = o(a), yielding 5(0) = and continuity of 



Notice that differentiability in 5 for F is assumed at the origin only, 
and that no regularity is asserted for 5(-) except at the origin. Assum- 
ing further regularity in a, for instance Holder continuity, we would obtain 
Holder continuity at a = of <$(•) by the same proof, incorporating the new 
bound |?t.2 (5, a) | < C|a| a coming from Holder continuity to work on ball 
B(0,2C\a\ a ). 

2.3.1 Application to bifurcation 

Lemma 12.271 is precisely what is needed to apply the Lyapunov-Schmidt 
method in finite-parameter bifurcation analyses. Consider a continuous curve 
0, as in (|2.25|) . taking for example u)(e, a) = 0. 

Proposition 2.28. Suppose that Ni,N 2 satisfy f)2.3p . (12. 4p . that R satisfies 
( 12.51) . (12. 9p . and that Assumptions 12.41 and 12.241 hold. Suppose further that 
a G R. Then, there exists a map a \— > 60(a), £0(0) = 0, continuous at 0, 
defined on a neighborhood ofO, and such that 



in a neighborhood of the origin in M a , where B is the continuous map given 
by Proposition \2.9[ 



2 sup 1 712(5, a) I = o(H), 

|S|<c 




□ 



(/, 50Ma) 5 a, B(e s (a), a, 0(a))) = 0, 
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Proof. The proof of Proposition ^. 91 with the exception of its last paragraph, 
still holds. Substituting continuity for Lipschitz continuity in the argument 
of that final paragraph, we obtain a continuous (in £>i-norm) map (e, a) i— > 
B := B(e, a, 0(a)) satisfying (j2.17j) . and such that g(e, a, B) = 0, for |e| + \a\ 
small. Then, following the notations of the proof of Proposition l2.12l consider 
the continuous map (e, a) i— > f(e, a), /(0, 0) = 0. By assumption on R, f(e, 0) 
is differentiable at e — 0, and d e f(0, 0) = 1. We can apply Lemma 12. 2 7\ and 
a solution to / = yields a solution to / = 0. □ 

3 Linearized bounds 

In the remainder of the paper, we apply the framework of Section [2] to the 
problem of bifurcation of viscous shock waves, first in the one- dimensional 
and then in the multi-dimensional case. We begin by assembling the relevant 
linearized estimates on the one- dimensional flow near a viscous shock. 

3.1 Hypotheses and notation 

As described in Section ll.4[ consider a one-parameter family of standing 
viscous shock solutions 

(3.1) u(x, t) = u £ (x), lim u £ (z) = u± (constant for fixed e), 

z^±oo 

of a smoothly- varying family of conservation laws 

(3.2) ut = T{e, u) := u xx - F(e, u) x , u E R n , 
with associated linearized operators 

(3.3) L{e) := dTjdu\ u= ^ = -d x A £ (x) + 8%, 

A £ (x) := F u (e ) u £ (x)) ) denoting A £ ± := lim 2 ^ ±00 A e (z) = F u {e,u e ±). Profiles 
u £ satisfy the standing-wave ODE 

(3.4) u' = F(s,u) -F(e,u e _). 

We recall the standard assumptions of [ZH, Z4, ZlJ: 
(HO) F E C k , k > 2. 
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(HI) o~(A e ± ) real, distinct, and nonzero. 

(H2) Considered as connecting orbits of (|3.4p . u e are transverse and 
unique up to translation, with dimensions of the stable subpace S(A e + ) and 
the unstable subspace U(A e _) summing for each e to n + 1. 

Remark 3.1. Condition (HI) implies that u £ ± are nonhyperbolic rest points 
of ODE ([32D, Ma(F u (e,u £ ±)) ^ 0, whence, by standard ODE theory, 

(3.5) \d£(u e -u e ± )(x)\<Ce- r >W, 0<£<k + l, 

for x ^ 0, some rj, C > 0; in particular, |(M e )'(x)| < Ce^L Condition (H2) 
implies in part that u e is of standard Lax type, i.e., the hyperbolic convection 
matrices A e ± at ±oo have, respectively, p — 1 negative and n — p positive 
eigenvalues for 1 < p < n, where p is the characteristic family associated 
with the shock. For further discussion, see |ZH| |ZT| [HZJ. 

To (H0)-(H2) we adjoin the generalized spectral (i.e., Evans function) 
condition: 

(T> e ) On a neighborhood of {9?A > 0} \ {0}, the only zeroes of D are 
(i) a zero of multiplicity one at A = 0, and (ii) a crossing conjugate pair of 
zeroes A±(e) = 7(e) + ir(e) with 7(0) = 0, <9 £ 7(0) > 0, and r(0) ^ 0. 

Lemma 3.2. Conditions (H0)-(H2) and (T> e ) are equivalent to conditions 
(P)(i)-(iii) of the introduction (stated in Section U..2.3}) together with F G C k , 
k > 2, simplicity and nonvanishing of a(A±), and the Lax condition 

(3.6) dim S(A £ + ) + dim U(A £ _) = n + 1, 

with u e (linearly and nonlinearly) stable for e < and unstable for e > 0. 

Proof. By (HO), (H2) and standard ODE theory, solution u e of is C 3 
in x, C 2 in (x,e), and decays at exponential rate in first two derivatives to 
endstates u e ± as x — > ±00. Thus, L(s) is asymptotically constant-coefficient, 
and its essential spectrum may be computed as in (jl.8D - (jl.9p . Section [1~2| 
to lie entirely within a parabolic region: 

a css (L(e)) C {3?A < -^|SA| 2 }, > 0, 

verifying (P)(i). On the other hand, computation (11.91) shows that (P)(i) 
implies (H2) in the generic situation that cr{A±) is simple and nonvanishing. 



45 



Likewise, recalling [AGJl [GZJ that to the right of the essential spectrum 
boundary, zeroes of D £ correspond to eigenvalues of L(e), we find that (V £ )(i) 
and (ii) are exactly (P) (ii) and (iii). Finally, recalling the stability criterion 
(see, e.g., [ZHl IZ4l IMaZ3l IMaZ4j ) that, under (H0)-(H2), linear and nonlinear 
stability of u £ is equivalent to 

(Z>) On {5iA > 0}, D £ has a single zero at A = 0, 

we find that u £ is stable precisely for e < 0. □ 

Remark 3.3. The undercompressive shocks studied in |HZ] give an example 
of shocks satisfying (P) but not (13. 6p . Their weaker decay properties are 
insufficient for the convergence argument of Section 11.6.41 see Remark 13.101 

Remark 3.4. Under (H0)-(H2), condition T> £ (i) is equivalent to (£ £ , (u £ )') = 
£ £ -{u%-u £ _) ^ for £ £ (constant) orthogonal to S(A £ _)UU(A%) [GZJ. Under 
the normalization (£ £ , (u 6 )') = 1, operator U £ Q f := (u £ )'(£ £ , f) plays the role 
of a "generalized spectral projection" onto KeiL(e) = Span{(tt e )'}, and £ £ 
the role of a generalized left eigenfunction [ZHj . Note that £ £ lies outside the 
domain of 1T5, a consequence of the absence of a spectral gap. 

Finally, we introduce the Banach spaces B\ = L 2 , £> 2 = c^L 1 fl L 2 , 
X l = {f: \f(x)\ < C(l + Ixl)- 1 }, 

and 

x 2 = d x {f: \f(x)\<C(i + \x\)- 2 }nx 1 , 

equipped with norms = ||/|| L 2, \\d x f\\B 2 = ll/IU 1 + ll^/IU 2 > 

\\f\\ Xl = ||(1 + \x\)f\\ L oo, and ||^/|U 2 = ||(1 + \x\) 2 f\\ L ~ + \\d x f\\ Xl , 

where d x is taken in the sense of distributions. By inspection, we have that 
£> 2 C £>i, X 2 C Xi, Xi C £>i, X 2 C £> 2 , and the closed unit ball in Xi is 
closed as a subset of B\. 

3.2 Projector bounds 

Lemma 3.5. Under assumptions (H0)-(H2), {T> £ ), associated with eigenval- 
ues A±(e) of L(e) are right and left eigenf unctions (f) £ ± and $ £ ± e C k (x,e), 
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k > 2 as in Assumption (HO) from Section [XT], exponentially decaying in up 
to k derivatives as x — > ±00 , and L(e) -invariant projection 

(3-7) Wf:=Y,^Wpf) 

i=± 

onto the total (oscillatory) eigenspace S e := Span{0±} ; bounded from L q or 
Ei to W 2 ' p fl X 2 for any 1 < q,p < oo. Moreover, 

(3.8) </4 = ^.$ £ ± , 

with $ e G C* +1 exponentially decaying in up to k + l derivatives as x — > ±oo. 

Proof. From simplicity of A±, we obtain either by standard spectral pertur- 
bation theory [KatJ or by direct Evans-function calculations |GJ1| IGJ2| IZH] 
that there exist A±(-), <f)±(-) € L 2 with the same smoothness C 2 (e) assumed 
on F. The exponential decay properties in x then follow by standard asymp- 
totic ODE theory; see, e.g., [GZ[ IZlj . Finally, recall the observation of [ZHJ 
that, by divergence form of L(e), we may integrate L(s)<j) = A</> from x = — oo 
to x — +oo to obtain A f^°° (p(x)dx = 0, and thereby (since \± ^ by as- 
sumption) 



(3.9) / (p±(x)dx = 0, 

J — oo 

from which we obtain by integration (13.81) with the stated properties of $±. 
From f)3.8p and representation (13. 7p . we obtain by Holder's inequality the 
stated bounds on projection n e . □ 

Defining IF := Id — n e , S e := RangelF, and L(e) := L(e)Il £ , denote by 

(3.10) G(x,t;y):=e Lie)t 5 y {x) 

the Green kernel associated with the linearized solution operator e Lt of the 
linearized evolution equations Ut = L(e)u, and 

(3.11) G(x,t;y) := e L{£)t fl £ 5 y (x) 

the Green kernel associated with the transverse linearized solution operator 
e £(e)*n £ . By direct computation, G = O + G, where 

(3.12) 0(x,t;y) := e W£)+ir(e)) V + (x)0' + (?/) + e { ^ £) ~ iT{£))t <p.{x)^_{y). 
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3.3 Short time estimates 



Lemma 3.6. Under assumptions (H0)-(H2), (T> £ ), for < t < T , any fixed 
T > 0, \e\ < So, there exists C > 0, depending only on Eq and T, such that 



(3.13) 
(3.14) 

and 

(3.15) 
(3.16) 

Besides, 

(3.17) 
(3.18) 



\\e L ^Wf\\ Bl < C\\f\\ Bl , 

\\e L ^dj\\ B2 < cdl/H^+r^H/llBj, 

\\e L ^Wf\\ Xl < C\\f\\ Xl , 
le^WdJU, < Ct-^ 2 snp(l + \x\ 2 )\f(x) 



||d £ (e Z ^fF)/|| 6l < C||/|| Sl , 
\d £ (e L ^) dx f\\ B2 < C(\\f\\ L . + \\f\\ Bl ] 



Proof. From standard parabolic semigroup bounds 

\e L(£)t \ L ^ L2 < Ct, \e L ^d x \ L 2^ L 2 < Cr 1 ' 2 , 

and properties 

(3.19) e L(£)t = e L(£) 'n £ + e L{e)t fl £ , 
and Ull 5 ^/!!^ < \f\ L 2, we obtain 

(3.20) \\e L ^d x f\\ L 2, || e z ^n e aj|| L2 < cr^WfU*, 

and, in particular, (13.1 3p . Likewise, we may obtain integrated bounds 



(3.21) 



Me)t d x f 



L 1 



e L{e)t fl E d x f 



L 1 



<c\\fl 



using the divergence form of L(e), by integrating the linearized equations 
with respect to x to obtain linearized equations U t = C(e)U for integrated 
variable 



U(e 1 x ) t) :-- 



u(e,z,t)dz, u(-,t) := e L{e)t d x f, 
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with linearized operator C(e) := — A e {x)d x + d 2 of the same parabolic form as 
L(e), then applying standard parabolic semigroup estimates (alternatively, 
pointwise Green function bounds as in Proposition 13. 7p to bound 

\\U(;t,e)\\ L i = \\e c ^f\\ Ll <C\\f\\ L i, 

the e Lt Ii £ bound then following by relation e L ^ 1 = e L ^ £S>t U £ +e L ^ t Il £ together 
with \\JWd x f\\ L i < \ f\ L i. Combining (15301) and (I52B . we obtain flBTHj) . 
Bounds f)3.15p and (I3.16P follow by identical arguments, substituting for L 2 
and L 1 the weighted spaces ||(1 + \x\) ■ and ||(1 + \x\) 2 ■ 

Bound (13.171) follows by the observation that v := d £ e L ^ t f solves vari- 
ational equation v t — Lv = (d £ L)e Lt f with initial data v t =o = 0, hence, by 
Duhamel's principle, 

v (t) = [ e Lit - s) {d £ L)e Ls fds. 
Jo 

By ( 13.131) . therefore, noting that \\(d £ L)g\\ Bl < C(\\g\\ Bl + \\d x g\\ Bl ), 



(3.22) \\v(t)\\ Bl <C Al + S - 1/2 )||/|| Bl ^<C 2 

Jo 



Bl- 



and (I3.17P follows from (I3.22p . the ^-regularity of IP stated in Lemma |3.5[ 
(13~T3|) and fl3~T9|) . 
Finally, let 

W(e, x, t) := I w(e, z, t)dz, w(; t) : = {d £ e L ^ l )d x f. 

J — oo 

Using again the divergence form of L, we find that W satisfies a parabolic 
equation 

d t W-CW= [ (d £ L)e tL d x fdx = d £ (d u F(e,u £ ))e tL d x f J 

J — oo 

where £ is of the same form as L; this implies 

\\W(t)\\ L i < [ \\d e (d u F(e,u*))\\ L oo\\e aL d x f\\ L id3, 



from which we obtain 
(3.23) \\W(t)\\ Ll <(C I s~^ds)\\f\\ L i, 
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Here, we have used the standard bound ||e* L 9 x /||LP < Ct _1 / 2 ||/||iP, 1 < p < 
oo, following from pointwise (high-frequency) resolvent estimates obtained 
by asymptotic ODE theory [Sj IPWj IZH] . To establish (13.181) . it remains to 
bound ||u>||l 2 - This is done as in the proof of (13.1 7ft : 



\w(t)\\ 



LP 



\d e L)e Ls d x fds 



LP 



[3.24) 



^ e L ^d x (d £ (d u F(e, u e ))e Ls d x f) ds 

<C f\t-s)- l l 2 \\e Ls d x f\\ Lv ds 
Jo 

{t _ s) -l/2 s -l/2 ds y ifhp 



Li' 



< c 



Applying to p = 2 gives the result. Estimate (13.181) follows from (I3.23p . 
(13^4j) and fl339|) . □ 



3.4 Pointwise Green function bounds 

Supressing the parameter e, denote by af, if, and rj the eigenvalues and 
left and right eigenvectors of A± = F u (e,u £ ± ). 

In the following Proposition, t] and M denote positive constants, 

\v\ 



\ a k 




(3.25) z%(y } t):=af[t 
and 

(3.26) P%(x,t;y):= 

are approximate scattered characteristic paths and effective diffusion rates 
along them, scattering coefficients [<4' l _], i — ~ , 0, + are constant, x ± denotes 
the positive/negative part of x, X{t>i} is a smooth cutoff function in t, iden- 
tically one for t > 1 and identically zero for t < 1/2, and indicator function 
X{| "t|>|y|} i s one f° r \ a k^\ — \v\ an d zero otherwise. 

Proposition 3.7 ( [ZHl [Z4l [MaZ3] ). Under assumptions (H0)-(H2), (V E ), 

(3.27) G = £ + S + K, 
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G as in (13.111) . where, for y < and all t > 0: 
(3.28) 

£(x,t;y) := £ [<_]£/' W (errfn (^^) " errfn (l^S* 

a->0 

and 

S(x,t; y) 

■= X{t>l } E r fc -/ fc - < (4 7 rt)- 1 /2 e -(-,- fc -*) 2 /4. 

a~<0 

+ E r 4 -t , (*rt)-^ e -«^« , '« (-^j) 

(3.29) a fc>° 

a fe >0, <0 

+x,.m, E Ki]r;/ t -'(4^jr 1 ' 2 e - (i -<« |2/ <> , (-^:) 



a k >0,aj>0 



denote excited and scattering terms and 
(3.30) 

TZ(x,t;y) = 0(e-^ x ' yl+t) ) 



n 

+ I] ((* + l)- 1/2 e-^ + + e-" 1 " 1 ) t" 1/2 e- 
fc=i 

+ E ^KT*|>|,|}0((t + l)- 1/2 t- 1/2 )e- 



(x-y-a~t) 2 /Mt 



(x-z- k ) 2 /Mt e - V x+ _ 



a fe >0, o i <0 



+ E x { |. l -«i, l „| ) o(( i +i)- i /v^ )e -(- 



a k >0,aj>0 
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n 

K t (x, f ; y) = 0{e-^ x ~ y \ +t) ) + E ° ((* + l)~ 1/2 e~' ,s+ + e ~ T/W ) 



k=i 

\2 , 



x (t + l)V2 t -3/2 e -(^- ?/ -a fc t) 2 /Mt 

(3-31) + £ X {|afc -„ l?/l} 0(t- 3 / 2 )e-^^) 2 /-e-- + , 



a fc >0,a J <0 



4- V v n(t- 3 P)p-( x - z tk) 2 / Mt p- r > x ~ 



a,;>0,a+>0 



(3.32) 

K y (x, t; y) = d t r(x, t; y) + 0( e -tfl*-wl+*)) 

n 

+ E + l)" 1/2 e-" x+ + e-" |a;| ) r i e -(-^*) 2 /^ 
fc=i 

Qfe >0,oJ<0 



a fe >0,aj>0 



r(x,*;j/) = 0(e-" ls/| (t + I)" 1 / 2 ) 
(3.33) x ( E e-^-^ 2 ^* + e -(*-w-«£*) 2 /A«^ . 



a k <0 



a t %(x, t; y) = 0(e-^ x - yl+t) ) + J2°(( t + ^)~ 1/2 e~ vx+ + e~ r ' lx ^ 

k=i 

x (t+ if/h'^+^^e'^-y-^^ 2 / Mt 

+ E X { | afc -,|>|,| } 0(t- (fc+3)/2 )e- ( ^ z - )2/Mt e-- + , 
(3.34) ^ >0 '^ <0 



4- V V f)|Y-(fc+3)/2\ p -(x-z+) 2 /Mt -r,x- 

+ 1^ X {\Ht\>\y\} U \ t ) e e 



a k >0,aj>0 



+ 0(e^l(t + l)-( fe + 2 )/ 2 ) 
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k > 1, a faster decaying residual. 

Moreover, for y < 0, k > 1, some M > 0, 
(3.35) 

n 
k=l 



x (t + l^/2 t -(k+l/2) e ~(x-y-a k t) 2 /Mt 

+ E X {|afc - t| > |j/|} 0((t + l)-^ 2 )e-^) 2 ^e-- + 



a fc >0, <0 



+ E X(K.|>l«0(«+l)- ( ' +2,/2 ) 



(fc+2)/2^ e -( x -z+) 2 /Mt e - V x- 



a k >0, aj>0 



+ 0(e-^l(t+l)-( fc+2 )/ 2 ) 

x ^ ^ e -(x-!,-«^) 2 /Mf + e _ (a: _ 2/ _ a + t )2 / ^y 

a~<0 a fe>° 

Symmetric bounds hold for y > 0. 
Proof. Evidently, it is equivalent to establish decomposition 

g= o+s+s+n 

of the full Green function. This problem has been treated in [ZH| IMaZ3j. 
starting with Inverse Laplace Transform representation 

(3.36) G(x, t] y) = e Lt 5 y (x) = j e A '(A - L(e))~ 1 5 y (x)d\ , 

where 

r:=<9{A:SRA<77i-?7 2 |3fA|} 

is an appropriate sectorial contour, rji, r] 2 > 0; estimating the resolvent kernel 
G\(x,y) := (A — L(e))~ 1 5 y (x) using Taylor expansion in A, asymptotic ODE 
techniques in x, y, and judicious decomposition into various scattering, ex- 
cited, and residual modes; then, finally, estimating the contribution of various 
modes to (I3.36P by Riemann saddlepoint (Stationary Phase) method, moving 
contour T to a optimal, "minimax" positions for each mode, depending on 
the values of (x, y, t). 



53 



In the present case, we may first move T to a contour T' enclosing (to the 
left) all spectra of L(e) except for the crossing pair A±(e), to obtain 

G(x,t;y) = I e xt (\- L(e))- 1 d\ + V Residue A (£) (e A '(A - L(£))~\(:r)) , 
Jr> j=± 

where Residue i.( e ) (e A< (A — L(e))~ 1 5 y (x)) = 0(x,t;y), then estimate the re- 
maining term <£ r , e A< (A — L(e))~ 1 dX on minimax contours as just described. 
See the proof of Proposition 7.1, [MaZ3], for a detailed discussion of minimax 
estimates £ + S + 71 and of Proposition 7.7, [MaZ3j . for a complementary 
discussion of residues incurred at eigenvalues in {3? A > 0} \ {0}. 

We have repaired in (13.301) . f!3.32|) two minor omissions in the bounds 
of |MaZ3j . pointed out already in [R], IHR[ IHRZj . Specifically, (i) in the 
first term on the righthand side of (13.301) . (13.321) . we have replaced the 
term 0(e~ vt e~\ x ~ y \ / M *)) appearing in |MaZ3l [ZTj . with a corrected version 
0(e-^\ x -y\ + V), and (ii) we have added a missing term 
(3.37) 

d t r = d t O(e-^(t + I)- 1 ' 2 ) ( e- {x - y -^ t?IMt + £ e -i*-v-°it) 2 /m^ 

a~<0 a fe >0 

which may alternatively be estimated as in [Rj IHRl IHRZ] as 

(3.38) d t r = 0(e-^\t + l)- 1 ) ( ^ e -(-y-^r/Mt + ^ ^-y-a+tf /Mt^ 

a~<0 a t >0 

to the righthand side of (I3.32p . 

As discussed in [R], the first correction concerns only bookkeeping, and 
comes from the fact that approximate Gaussians appearing in S decay only as 
e -mx-y\ j n f ar g e y anc | no ^. ag does total solution as 0{e~ r!t e^ x ^ / Mt )). 

This is only an artifact of our way of displaying the solution. The second 
is more fundamental, and deserves further comment. Namely, different from 
the constant-coefficient case, the (approximate) Gaussians appearing in the 
Green function for the variable coefficient case involve convection, diffusion, 
and also directional parameters that depend on x and y through u e (x). Thus, 
x- and ^/-derivatives do not always bring down additional additional factors 
i -1 / 2 of temporal decay as for constant-coefficient Gaussians, but may in- 
stead, through derivatives falling on spatially-dependent parameters, bring 
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down factors e - ^' or e - ^' of spatial decay; see [ZHj [Z4] for further discus- 
sion. The neglected term (13.381) . which is harmless for the stability analyses 
of [Rl IMaZ2l IMaZ4l IZlj , may be recognized as a term of this form. 

Unfortunately, term (I3.38p . though harmless in the stability analysis, is 
unacceptable for our bifurcation analysis (see Remark [3. lOj) . and thus must 
be treated in a different way. Reviewing the analysis of [MaZ3j . Proposition 
7.7, we see that term (I3.38P results from a term in (I3.36P of form 

(3.39) d y j e xt \<p(x)?p(y)d\, 

where Ti is a fixed finite arc of a larger contour T, ip(y)i d y ip(y) ~ e _r? ' y ' 
is a fast-decaying mode and 4>{x) ~ e ci3?A+c 2 5ft(A 2 ) j g a giow-decaying mode; 

the factor A reflects the fact that this is a Taylor remainder term in the 
low-frequency expansion about A = 0. The "old" estimate (I3.38P makes use 
of the fact that each factor of A introduces (by a simple scaling argument) 
an additional factor {t + l) -1 / 2 of temporal decay in the Riemann saddle- 
point /Stationary Phase estimate from which the final bounds are obtained. 

The "new" estimate (I3.37P follows by a still simpler argument, making use 
of the fact that a factor of A corresponds identically to time-differentiation: 

d y (f e xt \<p(x)ip(y) d\ = d t (f e xt d y <j){x)${y) d\, 

and estimating j> Fi e M d y 4>(x)ip(y) d\ as before. 

Conversely, a time-derivative on the Green function can always be traded 
for a factor of A in the integrand of f!3.36p . yielding faster temporal decay by 
factor (i: + l) x 'H~ l ~ t^ 1 ^ 2 . We use this observation to generate estimates 
f l3T3T|) . (I3734D . and (ET35]) not stated in [MaZ3llZl] . Note that conglomerate 
high-derivative bound (I3.35j) does not contain the "bookkeeping" error term 
Q( e -v(\%-y\+t)^ deriving from our way of decomposing the solution, and is in 
fact much easier to get than the more detailed individual bounds, requir- 
ing only modulus-based stationary-phase estimates like those used to bound 
residual 1Z. □ 

Remark 3.8. The effect of t-derivatives for variable-coefficient equations 
in bringing down additional factors of temporal decay is in contrast to that 
of x- and ^/-derivatives, which eventually bring down also e _r? ' z ' and e~ Tl ^ 
factors not adding extra decay [ZH4 IHZj . 
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Remark 3.9. Integrating the pointwise bounds of Proposition ^. 71 we obtain 

|| (G - £)y\\ L 2(x) < Ct _3/4 , || (G - £)y t \\L 2 (x) < CiT^ 4 , 

similarly as for the convected heat kernel K in the analysis of Section 11.6.41 

Remark 3.10. Note that G bounds indeed agree with the model given 
in Section 11.2.11 Specifically, £ terms are exactly superpositions of terms 
of form J. Likewise, S terms near their centers are well-approximated by 
moving Gaussians K, so satisfy similar bounds; see [Rj. Thus, the model 
analysis well- represents the primary terms £ and S in the transverse Green 
function expansion for Lax-type viscous shocks. 
Undercompressive shocks give a different model 

G ~ K(x, t; y)p(y) + J (x, t; y)q(y) 

with \p y \, \q y \ ~ e _r? ' y ', and thus 

(Kp) y ~ Ke- vlyl , 

insufficient for the argument. For further discussion, see |Z41 HZJ. A related, 
more subtle point arising also in the Lax shock case is that terms ~ e~ v ^\K y \ 
appearing in R y bound (13.381) are also not acceptable, since the e~ v ^ term 
is no help in convolution against an L 1 source, and \K y \ does not have the 
sign cancellation of K y . Our replacement with bound (13.371) amounts to 
showing that, even though the ^/-derivative no longer introduces cancellation 
for these terms, they already possess an additional t-derivative that leads to 
cancellation by another route. 

3.4.1 Parameter-dependent bounds 

Bounds involving e-derivatives do not appear in [ZH[ IMaZ3j . where depen- 
dence on parameters was not considered. Such bounds are necessary here 
in the uniqueness analysis. Recall, Lipschitz continuity with respect to the 
bifurcation parameter is needed in Section [2] to obtain an uniqueness result. 

Proposition 3.11. Under assumptions (H0)-(H2), (V e ), the residual 1Z, 



56 



introduced in Proposition 13.71 satisfies 
(3.40) 

d e Ky(x, t; y) = d t d e r(x, t; y) + d t r(x, t; y) + 0(e~^ x ~y^) 



fc=i 

+ E X {la , tl> _ lyl} 0((t + l)-^)e 

a k >0, o^. <0 
a fe >0,a+>0 



-l^-ix-y-a-t^/Mt 



-(x-z jk ) 2 /Mt e -r\x^ 



(x-zf k ) 2 /Mt e -rix- 



(3.41) 



(3.42) 



f(ar, t; y) = E ((* + l) 17 ^ - *** + (t + l^""'*') r i e -(*-»-« fc *) a /^t 



fc=i 



+ E Xfla^l^OCCt+l) 1 ^" 1 )^ 



(a=-^fc) 2 /^* e -^ + 



a fc >0, <0 



a h >0,aj>0 



57 



Moreover, for y < 0, k > 1, some M > 0, 
(3.43) 

n 

dfd e G y (x,t;y) = Y,° (r 1/2 e^ + +e-^l) 



k=l 



x (t + iY k+i y 2 t^ (k+1) e~ {x ~ v ~ a ^ t)2/Mt 

+ E X fl . :t |>l»| ) 0((« + l)- ( ' +1,/2 )^ ( "-^'' /M, e-- + , 



(fc+l)/2\ -( x - Z +f/Mt- V x- 



+ E x (Kl| > lall o((t + i)-"= + «/ 2) 

a fc >0, a+>0 

+ 0(e-^l(t + l)-( fc+2 )/ 2 ) 

x ^ ^ e -{x-y-a^t) 2 /Mt + ^ g-^-a+t) 2 /^ 



a+>0 



Symmetric bounds hold for y > 0. 



Proof. The above bounds may be obtained in similar fashion to those of 
Proposition 13.71 but require some additional discussion. Specifically, re- 
call in the construction of [MaZ3j that, in the crucial low-frequency regime, 
d y 1Z\(x, y) is expanded about A = as a sum of terms of the form 



(3.44) d y j e M d(\)(j)(x,\)ip(y,\)d\, 

d analytic, where 0, ip are of form 

0=(/ + ± (A,x))e^N(A), 

(3.45) . ... 

j> = (I + Q ± {\,y))<* x »w{\), 

with //, v, v, w analytic, and 

KdlQ ± (x,X)\<Ce-^, x^O 
\d{d y e ± (y,\)\<Ce-^, y^O 

for < i < k, < j. These are then estimated by the Riemann Saddle- 
point method as described above. Contributions in mid- and high-frequency 
regimes are of negligible order 0(e _r? (l a; l+l2/]+*)) and can be handled by simpler 
estimates [ZHllMaZ3j . 
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To treat the parameter-dependent case, we have only to observe that 
fl3.44p - fl3.45p still hold, but with d, /i, u, v, w now depending on e as well 
as A, and 0± and 0± depending on e as well as (x, A) and (y, A), where the 
regularity inherited from (HO) is C k in both x and e, k > 2 as in Assumption 
(HO) from the introduction of Section 13.11 and 9±, 9± by the same fixed- 
point/contraction mapping construction as in [MaZ3j satisfy 



(3.47) 



\d r x did{Q±(x,\,e)\ < Ce""'*', x^O 
\d r x dtd{e ± (y,X,s)\<Ce-^, y^O 



for < i,j < k, < r. Further, just as in the parameter-independent case, 
growth rates /i, v are either fast modes > c > (resp. \u\ > c > 0) for all 
e, or slow modes \i = Xfio (resp. v = Xuq), where /io (resp. z/ ) is analytic in A 
and C k in e. Moreover, either e^ x and e vy are both decaying, ^ftfix, tflvy < 
for 9ftA > 0, or else fi = —u, with e^ x ~ y ^ decaying, ^Rfx(x-y) < for 9ftA > 0. 

We treat e-derivative terms case-by-case, depending on whether /i and 
v are fast or slow, and which factor the ^-derivative falls on in the Leibnitz 
expansion. If d £ falls on ®± or 0±, for example, then by (13.471) the result is no 
worse than in the undifferentiated case, and so these terms may be subsumed 
in the previously-obtained bounds (I3.32p for !Z y . Likewise, derivatives falling 
on v or w are harmless, yielding terms of form already estimated in (I3.32p . 
Derivatives falling on fast modes e^, e uy , or e^ - ^- necessarily decaying by 
the case structure described above- yield factors fiexe^, v s ye vy ', or fi £ (x — 
y ) e M*-v) bounded by terms Ce^/ 2 , CW 2 , or Ce^ x ~ y ^ 2 of the same order 
as the corresponding undifferentiated term. Therefore, these contributions, 
too, may be bounded by estimate (13.321) . 

It remains only to estimate terms for which the e-derivative falls on a slow 
mode e^, e uy , or e^ x ~ y ^- equivalently, e x>MoX , e Xu ° y , or e XfM °( x ~ y ^- yielding an 
additional factor of \d £ fiQX ~ Ax, \d £ i>Qy ~ \y, or \d E fio(x — y) ~ A(x — y). 
The A factor may either be accounted as a t-derivative, or else as an extra 
factor of (t + l) -1 / 2 decay arising through the stationary phase estimate. 
The remaining x, y, or (x — y) factors, when multiplied against the moving 
Gaussian bounds obtained by stationary phase, give a contribution of order 
(t+1). For example, (x—y) multiplied against e~( x ~ y ~ at ^ l Mt may be bounded 
by (t+1) times e~ ( - x ~ y ~ at ^ / 2M * ? since x — y ~ at+t 1 ^ 2 where the Gaussian is of 

significant size. Likewise, xory times the typical term X{\ a -t\>\ y \} e ~ ( ' X ~ Zjk ^ 
coming from the product of slow-decaying modes may be bounded by (i+ 1) 

times X{\ a -t\>\ y \} e ~~^ X ~ Zjk)2 /2Mt i since x ~ z jk + t 1/2 and \z jk \, \y/%\ < t. 
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Finally, y times the typical term e ^ y at ) 2 / Mt e coming from the product 
of a slow- decaying y-mode and a fast-decaying x-mode may be expanded as 

(y _ x \ e -(*-y-">t) 2 /Mt e -v\x\ _|_ xe -(x-y-at) 2 /Mt e -ri\x\ 

and each term estimated as above. Other cases go similarly 

Converting A to (t + l) -1 / 2 decay for terms belonging to d t r in (I3.32p . 
we obtain an estimate poorer by factor (t + l) _1 / 2 (t + 1) = (t + l) 1 / 2 for 
d £ r than for r itself. For all other terms of this type, we convert A into a 
time-derivative and group them into dtr, where f is bounded by (t + 1) times 
the remaining terms (besides d t r) in (13.321) . Combining these estimates with 
the previous ones, we obtain (I3.40p as claimed. 

Proceeding similary for the expansion of G, but converting all A factors 
to (1 + t) -1 / 2 decay, we obtain (I3.43p . similarly as in the proof of (I3.35p . □ 

Remark 3.12. The effect of e-derivatives, roughly, is growth by factor 
t 1 / 2 ; that is, differentiation with respect to e degrades our decay estimates. 
This makes derivative bounds much more delicate than sup-norm bounds to 
obtain. The t 1 / 2 growth rate may be understood by the formal computation, 
neglecting commutators, 

d £ e L ^ ~ t(d £ L)e L ^ ~ td x e L ^ ~ t l ^e L ^\ 

4 Return map construction 

We show in this Section how the question of existence and uniqueness of 
time-periodic solutions to (11.201) in a neighborhood of the stationary solution 
( 11.191) can be formulated in the abstract framework of Section [2l We use the 
notations set out in Section 13.11 in particular Assumptions (H0)-(H2) and 
(T> e ), and functional spaces X\, X 2 and B%, B2, and the short-time bounds of 
Proposition 13.61 

Given a family of stationary solutions u e (13. ip of the system (11.201) . and a 
family of dynamic solutions u e of the system (11.201) . define the perturbation 
variable 

(4.1) u(s, x, t) := u £ (x, t) — u £ (x), 
satisfying nonlinear perturbation equations 

(4.2) u t - L(e)u = Q(e,u) x , u(e, x, 0) = u (e,x), 
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where the nonlinear term 

(4.3) Q(e, u) := -F(e, u £ ) + F(e, u £ ) + F u {e, u £ )u £ , 
satisfies the pointwise bound 

(4.4) \Q{e,u)\ + \d £ Q{e,u)\<C\u\\ 

C being a nondecreasing function of || U fee . 
Decomposing 

(4.5) u = w + 4> £ + + W-(jf_ + v, 

where (with the notations of Section [372T) w + 4> £ + + W-4> £ _ = U £ u G £ 
IFw G £ e , and coordinatizing as (w,v), w := (w + ,wJ), we obtain 

w = (7(e) Id +t(e)J)w + N(e,w,v), 
(4.oJ ~ ~ 

v = L(e)v + H £ d x Q(e, w, v), 

*»;=(.', J), and 

iV(e, u>, u) := n £ 9 x .Q(e, + w_0L + d), 

«J, f ) := + w-<jf_ + u). 

From (14.31) . (I4.4p and the IF-bounds of Lemma [3.51 

(4.7) sup(l + |x| 2 ) |Q| < + I^IUJ 2 , \N\<C(\w\ + \\v\\ L 

X 

and, for 1 < q < oo, 

(4.8) \\d w Q\\ c{L ^ Lq) < C{\w\ + \\v\\ L ~), < C(H + \\v 
From (14.41) and Lemma I3.5[ for 1 < q < oo, 

(4.9) ||d e Q|| L9 <C<(M + IMI^) 2 - 

In (14.71) . (14.81) and (14. 9p . C is a non- decreasing function of \w\ + \\v\\l 
Now, truncate N, replacing it with 

(4.10) N(e, w, v) := N(e, w, v(w, v)), 
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where v(r,v) is the map 

(A „s ( a /^(Co , ^ )v(x,t), ifv(x,t)^0, 

(4.11) (x,t) i-> < |t>(:M)k» 

[ 0, ifv(z,t) = 0, 

Co being a positive constant and ip G M 1 being a C°° "truncation" function 
with 

f 1 z > 1, „ . f z > 1, 

(412) = ,i V 2, ^ = {1 2 i V 2, 

1^1 < 1) IV'I — 1- With these notations, \v\ < C\w\, and thus, by (14.71) with 

q = oo, 

(4.13) \N\ < C\w\ 2 , 

The map (w, v) eM 2 x [L q n L°°) (L 9 n L°°) is C°° at any (w, v) 

such that w^O, and infinitely Frechet different iable in v for any (w, v), with 
first partial derivatives 

d v v = i>(C ^)ld £{Lq) - C \w\i/(C l A^f^v, 
\v\ \v\ \v\ 6 

(with the convention (d v v ■ h)(x) := h(x) if v(x) = 0), and 

<9«,f = C — V> (Cof-^J, 

|w| \v\ \v\ 

(with the convention d w v(x) := if t> (x) = 0). In particular, 

(4-14) ||$>V|U(L«) < C, ||^|U(K2,L8) < C, 

for 1 < q < 00, uniformly in w, v, for w^O. 

Lemma 4.1. Under assumptions (H0)-(H2), (T> £ ), as a function from (w, v) e 
M. 2 x L 9 , /or any 1 < q < 00, N is C k , k > 2 as in Assumption (HO) from 
Section |3~TI /or w away from 0. At w = 0, N is C 1 . Moreover, 

(4.15) l^iVl^^) + ll^iVH^L,^) < C\w\. 
Besides, N is Lipschitz in e, with the bound, 

(4.16) \d £ N\ < C\w\ 2 . 

In (14.151) and (14.161) . C is a nondecreasing function of \w\ + 
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Proof. Estimate (I4.15P is immediate, using the chain rule, (14. 7p and (I4.14p . 
Estimate (14.161) follows from (14. 9 p and the e-regularity stated in Lemma 
1331 □ 

Consider the truncated system 

^ 417 ^ «; = (7(e)Id + r(e) J)w + N(e, w, v), 

v = L(e)v + f[ e d x Q(e, w, v), 

Proposition 4.2. Under assumptions (H0)-(H2), (T> e ), for < t < T, 
any fixed C\,T > 0, some C > 0, and \a\, ||6||xi> | £ | sufficiently small, 
system (I4.17P with initial data (w ,v ) = (a,b) possesses a unique solution 
(w,v)(e,a,b) e L°°([0,T],M 2 x X x ) that for a ^ is C k+1 in t and C k in 
(e,a,b), k > 2 as in Assumption (HO) from Section l3~Tj with respect to the 
weaker norm B\, and for a = is C 1 in t and Lipschitz in (e,a,b) with 
respect to E\, with 

C~ x \a\<\w{t)\<C\a\, 

(4 ' 18) IK*)lk<c(IHk + H 2 )- 

In particular, for ||fe||xi < Cil a l> all < t < T , 

(4.19) \Ht)\\x! < C\w{t)\. 
Besides, for all < t < T , 

(4.20) \d a w(t)\ + \\d a v(t)\\ Bl + Wdbwit)^^ + ||^(*)|U (Bl) < C, 
and 

(4.21) \d e w(t)\ + \\9ev(t)\\ Bl < C(\a\ + ||6|| Bl ). 

Proof. Existence and uniqueness follow by a standard Contraction-mapping 
argument, using Duhamel's Principle to express (14.171) as 

w(t) = W(t) := e Me)+Ti£)J)t a+ f e Ms)+T ^ )J){t - s) N(e, w, v)(s)ds, 

(4.22) Jo 

v (t) = V(t) :=e L{£)t b+ / e L{£Kt ~ s) fl £ d x Q(e,w,v)(s)ds. 

Jo 
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Consider indeed the map T e ^b : M 2 x Xx — > IR 2 x X 1; which maps (w,v) 
to (W, V), where W, V are defined in Estimates (jOl . (Q5II together 

with short-time bounds (I3.15p . (13.161) show that given any T > 0, for |e| < e , 
one can find p(eo,T) > such that, for \a\ + \\b\\x 1 < p, $e,a,b maps a closed 
ball in L°°([0,T],M 2 x X x ) to itself; the norm in L°°([0,T],R 2 x Xi) being 
defined as sup Q<t<T (\w(t)\ + Estimates (14. 8p . (I4.15P together with 

short-time bound (13.141) show that, provided that p is small enough, $ £ja ,b is 
moreover a contraction in that same ball of L°°([0, T], M 2 x Xi) with respect 
to the weaker norm sup 0<t<r (|w(t)| + whence we obtain existence 

and uniqueness by the bounded in strong norm/contractive in weak norm 
principle described in Section [2J 

Bound (14.181) (i) follows readily from (14.221) by Gronwall's inequality, using 
the "decoupled" estimate \N\ < C\w\ 2 for N. Bound (14.181) (ii) is a standard 
Gronwall bound that requires no comment. Combining (I4.18p (i)-(ii). we 
obtain evidently (14.191) for \a\ sufficiently small. 

Lipschitz regularity with respect to a, b of the fixed point of T e ^, still 
denoted (w, v), is a consequence of the uniform bound |<9 a W^(t)| < C and 
of the uniform estimate II^V"^))^^) < C (itself a direct consequence of 
(13~T3|) V Estimate flCTj) follows. 

Similarly, estimate (I4.2ip follows from the uniform estimates |0 e VF(t)| < 
C\a\ (a consequence of (14.161) ) and ||^ e T^(t) < C(\a\ + ||&||bi) ( a conse- 
quence of (I3~T7D . dSHBD and 

Finally, (14.181) (i) shows for a ^ that w(t) remains bounded from zero, 
whence the solution is C k , k > 2 in (e, a, b) and C k+l in t, by the correspond- 
ing properties of the righthand side of (14.221) away from w = 0. □ 

Remark 4.3. Since \\v\\ Xl controls \\v the bound (14.191) allows us to 
justify a posteriori the truncation step provided that we can later find a 
solution with < C|a| 2 , and \a\, \\b\\x 1 sufficiently small. That is, such 

a solution satisfies not only the truncated system (14.171) but also the original 
equations (14. 6p . 

Given (a,b) 6 M. 2 x Xi, such that a / 0, let (w,v) be the solution of 
(14TT7|) over [0,T] such that (w,v)(t = 0) = (a, b). Estimate (OSjl fi) shows 
that w never vanishes. Considering a double covering in polar coordinates of 
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the equation in the plane £ e , we coordinatize w as (r, 9) G M x T, obtaining 

r = 7(e)r + A^ r (£, r, 0, t>), 

(4.23) 9 = r(e) + N e (e,r,9,v), 

v = L(e)v + Wd x Q(e, re ie , re~ ie , v), 

where 

N r := \(e- w 4>% + e* e 0L, N), N e := ^(e^ft ~ N). 
Observing, by (Q5jl (ii). (|Q5jl . (OS]l (i) and (P) (iii) , that 

(4.24) e = T(e) + O(|o|)^0, 

we find that all values of 9 G [0, 27r] are taken on for T chosen sufficiently 
large; say, > 27r/r(e). Thus, in searching for periodic solutions, we may 
without loss of generality fix 0(0) = 0, or w(0) = (a,0), a G M 1 , thereby 
factoring out invariance under translation in t. 

With this choice of initial condition, we obtain after a brief calculation 

r(t) = e 7(e) *a + f e 7(e)( *~ s) N r (s,r,9,v)(s)ds, 
Jo 

(4.25) 9{t)=r{e)t+ [ N e (e,r,9,v)(s)ds, 

Jo 

v {t) = e L{£)t b + [ e L(£ ^ t -^fl £ d x Q(e,re i6 ,re- ie ,v)(s)ds 
Jo 

for all < t < T. Note that we have reduced the dimension of a from two 
to one; we will follow this convention for the rest of the analysis. 

Lemma 4.4. Under assumptions (H0)-(H2), {T> £ ), for \e\, \a\, \\b\\x 1 suf- 
ficiently small, a / G R 1 , consider the solution (w,v) of (14.171) on < 
t < 47r/r(0) issuing from the initial condition (w,v)(0) = ((a, 0),6), and let 
w = (re t6 , re~ l9 ). Then, there exists a unique smallest T = T(e,a,b) > 
such that 

(4.26) 9{e,a,b,T(e,a,b)) = 2ir. 

Moreover, the function (e, a, b) \— > T(e, a, b) is Lipschitz in E\ norm, with 

(4.27) T(e,0,0) = 2n/r(e). 
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Proof. The function 9 is well defined and solves (I4.25p (ii). Estimate (14. 181) 
shows for (a, b) = (0,0) that (r, v) = (0,0) for all t, hence 9{t) = r(e)t, 
giving (14.271) . Local existence, uniqueness, and regularity of T satisfying 
(I4.27P and (I4.26P then follows by the Implicit Function Theorem applied 
around (e, a, b, T) = (0, 0, 0, 2vr/r(0)), using the fact that fl424j) holds for all 
t < 47r/r(0); the function 9 inherits the regularity properties of w described 
in Proposition 14.21 and T inherits the regularity properties of 9. Finally, 
(I4.24p and continuity of r show that the minimal solution of (14.261) must lie 
near 27r/r(0), yielding global uniqueness as well. □ 

Substituting t = T(e,a,b) in (14.251) . we may express the Poincare re- 
turn map (e, a, b) — > (a,b) := (r, v)(e, a, b, T(e, a, b)) as a discrete dynamical 
system 

a = R(e, a, b)a + N\ (e, a, 6), 

(4.28) 

b = S(s, a, b)b + N 2 (e, a, 6), 

of the form (I2T21) studied in Section [2J with a, e, N x G R 1 and b G X 1 , 
N 2 G X 2 , where 

R(e,a,b) := e^ e)T ^ a ' h \ 

U 29) 

5(6, a, 6) : =e L ( £ ) T ^), 

are primary and transverse linearized solution operators for one time-step 
T(e,a,b) of continuous system (I4.17p . and 

l>T(e,a,b) 

iVi(e,a,6):= / e 7(e)(T(e ' a ' 6) - s) iV r (£, r, 9, v)(s)ds, 

N 2 (e,a,b):= / e i(£)(T(£ ' a ' fe) - s) n e a x Q(e, re lS , re^ e , v)(s)ds. 

Jo 

Evidently, small-amplitude periodic solutions of (I4.17P with period T 
close to T(0, 0, 0) = 27r/r(0) are equivalent to fixed points of the Poincare 
return map (equilibria of (I4.28P ). Moreover, on the wedge {H&Hx! < Ci|a|} 
(see (14.191) ). these are equivalent to small- amplitude periodic solutions of the 
original (untruncated) system (13.21) . by Remark [4.31 

Lemma 4.5. Under assumptions (H0)-(H2) and (T> £ ) (which are stated in 
Section I3TT1 and, by Lemma WHX are equivalent to the assumptions of Theorem 
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II. 3j) . R satisfies (j2.5p and (I2.8j) : R(0, 0, •) zs differentiable for alls sufficiently 
small and satisfies (12. 91) ,- Ni satisfies (12.31) and (12. 61) / iV 2 satisfies (12.41) and 

Proo/. From Lemma 031 i? satisfies (1231) and (jZ5]l : from (I4.27[) and (P) (iii) 
(see section 03]), we find that R(e, 0, 0) = e 27r7 ^/ T ^ is C fc , k > 2 for £ 
sufficiently small, satisfying (I2.9p . Likewise, the bounds on N% follow from 
bounds (g3HD (i), (14201) and (OS!) , and ODE bound 

| e 7(e)(T( e> o,6)-s)| < C 

on the primary linearized solution operator. Note that we find a control of 
|iVi| and \d EtCt ,bNi\ by \a\ 2 and \a\ respectively, stronger than (12. 3p . 

Bound (j£lD(i) follows from (13TT6D . (14\7D and (14~T8|) ; (l2~4]) fii) follows from 
(13~T4l) . (Ol) . (IQOP and (Qgp . and finally flU]) is a consequence of (fXTSj) . 
Q3ZEU) , (1QT1) . (HU) and dUBD- Note that the bounds for d £)Cl)b N 2 involve 
||5xQ(P)||bi, where Q is evaluated at e,r,9,v, (re id , re~ td , v ) being the so- 
lution of (14.171) given by Proposition I4.2| and T = T(e,a,b); by parabolic 
smoothing, \\d x Q(T)\\ Bl < C(l + T^) || 6 || Bi . ' □ 

Remark 4.6. Note that we have made use in Proposition 14.21 and Lemma 
14.51 of parabolic smoothing, reflected by integrability of singularity t -1 / 2 in 
the bound on e^^ILc^. Compare the analogous computations in the proof 
of Theorem 1.3, [TZlj on construction of a center manifold in the context of 
model problem I. 

At this point, we have reformulated model problem II in the abstract 
framework of Sections 12.11 and 12. 21 using only the conservative structure of 
the equations, the elementary semigroup estimates of Lemma 13. 6| and direct 
computation. It remains to verify Assumptions 12.41 and 12.61 by Proposition 
12.51 essentially linearized stability estimates on Y^jLo^i f° r which we shall 
require the detailed pointwise bounds (I3.27I) - (I3.26I) of Section 13.41 



5 Proof of Theorem 1.3 



We start from the perturbation equations (14.21) . where the linear operator 
L(e) is assumed to satisfy assumptions (H0)-(H2) and (V e ) (which are stated 
in Section 13.11 and are shown in Lemma 13.21 to be equivalent to the assump- 
tions of Theorem 11.31) . 
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Consider the Poincare return map system (14.281) associated with the trun- 
cated system (14.171) . derived from (14. 2 p by a proper choice of coordinatization 
and truncation. Periodic solutions of (I4.17P are in one-to-one correspondence 
with fixed points of (I4.28j) . Lemma I4~5l implies that the basic assumptions of 
Section [2TT1 are satisfied by the linear and nonlinear terms in (I4.28p . Thus we 
have reformulated the question of the existence and uniqueness of periodic 
solutions to the perturbation equations (14.21) in the abstract framework of 
Section [2j So far we used only the conservative structure of the equations, 
the elementary semigroup estimates of Lemma 13.61 and direct computation. 

We show in Section 15.21 that the assumptions of Proposition 12.281 are 
satisfied by systems (14.281) . thus yielding existence of periodic solutions to the 
perturbation equations (14. 2p . We show in Section I5"T31 that the assumptions 
of Corollary 12.171 are satisfied by system (I4.28p . thus yielding uniqueness of 
the periodic solution. Verification of these Assumptions require the detailed 
pointwise bounds of Sections 13.41 and 13.4.11 In a first step (Section 15. ip . 
we verify these assumptions for a model linearized operator described in the 
introduction. 

We use in this proof the notations (in particular, the definitions of the 
functional spaces Xi,X 2 and 81,82) set out in Section [37T1 

5.1 Key estimates 

We return in more detail to the model analysis of Section 11.6.41 taking 



where K(x,t;y) := t~ 1 / 2 e~^~ y ~ a ^ 2//4i models the scattering term S and 
J(x,t;y) := u'(x)eiTin((—y — at)/2t 1 ^ 2 ) (the function errfn is defined in 
( I1.12p ) the excited term £ in decomposition (13.271) . 

The aim of this Section is to check Assumptions 12.41 for the model trans- 
verse operator S defined by 



for some fixed T > 0, where G is given by (15. ip . so that the iterated transverse 
solution kernel takes form 



(5.1) 



G(%, t] y) = K(x, t;y) + J (x, t; y) 




S 3 5 y (x) = K(x,jT;y) + J(x,jT;y). 
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By Proposition ^. 51 to prove Assumption I2.4[ it is sufficient to prove that the 
sequence ^=o^ * s conditionally convergent in C(X 2 ,Bi) norm, and that 
the limit operator Y^jLo^ * s bounded in C(B2,Bi) and C(X2,Xi) norms. 
We assume noncharacteristicity, a < 0. Define the kernels 

N N 

K N := ^K{x,jT;y), J N := ^ J(x, jT;y), 

3=1 3=1 

and associated operators 

tC N f(x):= K N (x,jT;y)f(y)dy, J N f(x):= J N (x, jT;y)f(y) dy. 



5.1.1 Scattering estimate (K term) 

Approximating the kernel K N by integral Kn(x, y, T) := T~ l K(x, t; y) dt, 
define the "continuization error" kernel 9n '■= Kn — Kn, and associated op- 
erators 

/+oo r+oo 
K N (x,y,T)f(y)dy, N f(x):= 9 N (x,y,T)f(y)dy. 
-oo J —oo 

Lemma 5.1. There exists C > 0, such that, for all T > Tq, y 6 M, for 

a = 0, 1 and < (3 < 3, 

\\dtdyK(- } T,y)\\ L 2 < CT^ 1+2a+2 ^\ 
In particular, for T >T , 

/*+oo /*oo 

/ \\K yy (;t;y)\\ L2 dt+ / \\K yt (-, t; y)\\ L 2 dt < CT~ l '\ 
Jt Jt 

Proof. Direct computations, and use of \\x a e x2 ^ ct |jx? — 0(t a//2+1 / 4 ). □ 

Lemma 5.2. For all i,!/Gl, all T > T , there exists C(T) > such that 

(5.2) \K{x,T;y\<C{T){l + \x-y\)-\ 
There exists C > such that, for all x € R, T > T , 

(5.3) j {\K yy {x,t ]y )\,\K yt {x,t;y)\)dt,<C{l + \x-y\)-\ 
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Proof. The bound O follows from K(x,T;y) < C(T)e^ x -^ 2 / CT . Direct 
computations give the existence of C > 0, and c , c\ > 0, independent of 
x, y, t, such that, 

(5.4) \Kyy\ < cr 3/2 e- {x - y - at)2/cot , \Kyt\ < cr 3/2 e- {x - y - at)2/cit , 

In the light of (15.41) . to establish (15.31) . it is sufficient to prove the uniform 
bound 

/+oo 
\ x - y \t- 3 / 2 e-( x -y- at ^ ct dt<c. 

Observing that 

\ x - y \ t -3/2 e -(x-v-at)*/Ct < ( | x _ y _ at y t l/2^-l e -(x-y-at)VCt 

+ \ a \ t -V* e -(*-y~*tr/ct 

< e -(x-y-at) 2 /C 2 t 



we find in turn that it is sufficient to show 

(5.5) / t-We-^-v-^^dt^C, 



which follows by 

t -l/2 e -{*-V-<*)>/Ct dt 



T 



y-l/2 e -( X -y-at)yCt dt 



+ 

'\x—y—at\<t/C J\x-y-at\>t/C 

/+oo r+oo 
\ X -y\-V2 e -(*~y~at?/C\ X -y\ dt+ / t -l,2 e -t/C» dt 
-oo JO 

□ 

Lemma 5.3. For some C > 0, all x e M}, 

/+oo 
(l + |a:- 3 /|)- 1 (l + |y|)- 2 dl/<C7(l + |x|)- 1 . 
-oo 
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Proof. Let c(x,y) := (1 + \x — y\) 1 (1 + \y\) 2 . Dividing into commensurate 
and incommensurate parts, we have 



c(x,y)dy = / c(x, y) dy + / c(x,y)dy 

' J\x-y\<\x\/C J\x-y\>\x\/C 

px(l+l/C) r+oo 

<C 2 (1 + \x\)- 2 dy + C 2 (1 + \x\)-\l + \y\Y 2 dy 

Jx(l-1/C) J-oo 

< c 3 (i + l^i)- 1 . 

□ 

Lemma 5.4. The sequence JC^ is convergent in C(B 2 , £>i), uniformly inT G 
(T , oo). Its limit K,^ is Lipschitz inT G (T , oo) as an operator in C(B 2 , B\), 
and is bounded in C(X 2 ,X 1 ) norm, uniformly in T G (T , oo). 

Proof. That /C^ belongs to C(B 2 , Bi), for all iV, is a consequence of Lemma 
15.11 To prove convergence in C(B 2 , B\), we use the identity, 



K y = a-\K t -K yy ), 

which implies 

/•NT 

(5.6) aTd y k N = K(NT) - K{T) - \ K yy (s)ds, 

Jt 

hence 

/+oo 
d y (K N (x,y) - K N+P )f(y)dy 
-oo 

is bounded by 

C sup (\\K(x, NT; y)\\ L2{x) + \\K(x, (N + p)T; j/)^ 
y v 

+ / 11^(0;, t; y) \\ L ^ X ) dt) 

I NT 



L 2 (x) 



which, by Lemma 15.11 is in turn bounded by C(NT) 1 ^ 4 ||/||l 1 - This shows 
that /Cat(<9 x /) is Cauchy in £>i, uniformly in T > T , for d x f G B 2 . By the 
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uniform boundedness principle, the limit /Coo belongs to £(£>i,£> 2 ). By the 
Fundamental Theorem of Calculus, 

N f (j+l)T pT 

(5.7) 6 N = K(NT) + ~y2 / K t (s)dsdt, 

1 i JjT Jt 

hence 

/NT 
\K yt (x,y)\dt. 

This bound, together with the above Lemmas, implies that 0^ belongs to 
£(£>2, £>i), and 

/•+oo f(N+p)T 

\\(e N - e N+p )d x f\\ L > < / / |K t (^;t)|U/(y)|^t 



oo 



/+oo 
(\\K y (-,NT;y)\\ L2x + \\K y (-,(N+p)T;y\\ Ll )\f(y)\dy 
-oo 



< C(iVT)- 1/4 



L 1 



using Lemma [5.11 again. Hence &N{d x f) is Cauchy in £>i, uniformly in T > 
T , for all d x f G £> 2 . Combining, we have /Coo = ft^ + 6oo G C(Bi, £> 2 ). 

The kernel of /Coo is given by the limit N — > oo in the right-hand side of 
(15.61) . where both terms are convergent in the sense of £(82, B\). In particular, 
given d x f G X 2 , 
(5.9) 

/poo p 
K(x,T;y)f(y)dy- J J K yy (x,t;y)f(y)dydt, 

and, by Lemma lo"T2"| 

11^00(9*/) < Csup(l + \x\) f(l + \x- yD^imi dy, 

x J 



which, by Lemma [573| implies that /Coo G £(X 2 ,Xi), and that ||/Coo||£(x 2 ,Xi) 
is uniformly bounded in (T , 00). 

Similarly, using (15.71) . were both terms in the right-hand side were shown 
to converge as N — > 00, 



\Qoo(d x f)\\ Xl <Csup(l + |a;|) / \K yt (x,t;y)\\f(y)\dydt, 

X 

72 



and Lemmas 15.21 and 15.31 then imply that 0^ G C(X 2 ,Xi) as well, and 
a uniform bound on ||OooC0IU(x 2 ,Xi)- Combining, we have /Coo uniformly 
bounded in £(X 2 ,Xi) norm. 

To establish Lipschitz regularity of /Coo, we compute 

||«9 T (T/Coo)(^/)|| L2 < y(||^(x,T;y)|| Li + ||^(x,T;y)|| LS )|/(y)|^, 

so that, by Lemma loTT| /Coo is Lipschitz in (T , oo) as an operator in C(B 2 , £>i), 
with Lipschitz semi-norm bounded by CT _1 (1 + H/CooH^^^^). 

Lipschitz regularity of requires a bit more care, and an additional 
observation of general use. Namely, for a function like K y for which higher 
i-derivatives decay successively faster, the total truncation error 0oo can be 
partially evaluated as a time-independent function plus an arbitrarily rapidly 
converging integral in time, by using successively higher order numerical 
quadrature formulae. 

The second-order trapezoid rule, 

rN+l 

(/(l)/2 + /(2) + ■■■ + f(N) + f(N + l)/2) - J fit) dt 

= E( t f\st)f"(j + t)dtds+ f [\st + s-t)f"(j + t)dtds) 
■ =1 Jo Jo Jo Jl 



gives 



d y e N = (l/2)(K y (T) + K y ((N + l)T) 

N+l 



(5.10) 



+ / (-sr)T 2 K ytt (x,((j + T)T);y)drds 

j=1 K Jo Jo 

+ So Si ^ sr + s -^ T2 ) K y^ x ^^ + T ^y^ dTds )- 

In (15.101) . taking the limit N — > oo in £(£> 2 ,£>i), then differentiating with 
respect to T, we find that, given d x f G B%, 

d T Qoo{d x f)(x) = K yt (x,T;y)f(y)dy 

N+l i s 

+ E/ (/ (-ST)(2TK ytt + (j + r)T 2 K yttt )(x,((j + r)T);y)f(y)dr 
j =1 Jo v Jo 

+ J\st + s- r)((2TK ytt + (j + r)T 2 K yta )(x, ((j + r)T);y)f{y) dr)ds. 
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Lemma [5.11 then implies that 0oo is Lipschtiz in (T ,oo), as an operator in 
£(£> 2 ,£>i), with Lipschitz semi-norm bounded by CTq 1 ^. 

□ 

Remark 5.5. Similarly, we could use Simpson's rule to express 

dyOoofay) = (l/Q)K y (x,T;y) + (5/Q)K y (x,2T;y)+o(£ \K ytm {x,y)\dt) , 

obtaining convergence at rate f \\K ytm (x,y) \\u2( x )dt ~ / t~ n ^dt ~ t~ 7 ^, 
and so on. Together with the favorable t- derivative bounds noted in (13.351) . 
this allows us to neglect continuization error for any practical purpose of 
determining convergence or boundedness. 

5.1.2 Excited estimate (J term) 

Approximating the kernel Jn by integral Jn{x, y, T) := T~ l J(x, t; y) dt, 

define the continuization error kernel ipN(x,y) := Jn — Jn, and associated 
operators 



J N f{x) 



J N (x,y,T)f(y)dy, V N f(x) 



ip N (x,y,T)f(y) dy. 



Lemma 5.6. The sequence Jn is convergent in C(X 2 , Bi), uniformly inT G 
(T ,oo). Its limit J^ belongs to £(82, Bi), is Lipschitz in T G (T , 00) as 
an operator in £(B 2 ,Bi), on<i is bounded in C(X 2 ,Xi) norm, uniformly in 

Te(T Q , 00). 



Proof. By ( 13.51) . Remark 13.11 

\J y (x,t;y)\ = \c 2 u'(x)K(0,t;y)\ < Ce-^\K(0,t;y)\. 
Given d x f G X 2 , the crude bound \K\ < C then implies 

/NT p 
J J y (x;t,y)f(y)dydt ^ < CNT\f\ Lh 

and Jn G C(X 2 ,Bi). To prove convergence in C(X 2 ,Bi), we use \f(y)\ < 
\\d x f\\x 2 {l + \y\y 2 , to bound 



(N+p)T 



NT 



J y (x,t;y)f(y) dydt 
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by 

*(N+p)T 



j (l + \y\r 2 K(0,t;y)dt)\\d x f\\ X2 , 

then split (1 + \y\)~ 2 into (1 + |y|) _3 / 4 (l + |y|)~ 5 ^ 4 , where the second factor 
is integrable, to bound (15. lip by 



r (N+p)T 

(5.12) C( / \\(l + \y\)-V 4 K(0,t;y)\\ L ~dt)\\d x f\\ X2 . 

v J NT ' 

For t large enough, 

(5.13) SU P (1 + | y |)-3/4 t -l/2 e -(,-ai)V4t < C r^ 2 (l + |i|)" 3/4 , 

y 

so that for N large enough, (I5.12p is controlled by C(NT)~ l / A \\d x f\\x 2 i show- 
ing convergence of J7jv in C(X 2 , £>i), uniformly in T e (T , oo). 
The bound 

sup ((l + |x|) J j u'(x)(l + \y\)- 2 K(Q,t;y)dydtj < oo, 

is proved similarly and implies that the limit also belongs to £(X 2 ,X L ), 
and that \\Joo{T)\\c(x 2 ,x 1 ) is uniformly bounded for T e (To,oo). Note that 
it is a consequence of (15.51) in Lemma 15721 that 

r+ao 

(5.14) sup/ \\J y (x,t;y)\\ L 2 dt < oo. 

This implies in particular that j7oo also belongs to C(B 2 ,Bi), with norm 
controlled by (I5.14j) . Finally, 

MTJoo){d x f) = - f J y (x,T;y)f(y)dy, 

which implies that in £(£>2, Bi), Joo is Lipschitz with respect in T e (T , oo). 
The treatment of the continuization error \1/ jv and of its limit is straight- 
forward; see Remark 15.51 The results for Jn and j7~oo then follow by J7at = 
Jn + ^n- □ 
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5.2 Existence 



The nonlinear terms in (14.281) satisfy (I2.3P and (12. 4p . The linear operator R 
in (14.281) satisfies (12. 5p and (12.91) . In the (r, 6, v) coordinatization, the variable 
r is scalar. Thus, with the following Lemma, we can apply Proposition 12.281 
to obtain an existence result. 

Lemma 5.7. Under assumptions (H0)-(H2) ; (T> £ ), the linear operator S 
defined in (I4.29P satisfies Assumption \2A\ and \2.2A[ 

Proof. Following the model analysis of Section 15.11 approximate T,n : = 
Ylf=i & by its continuization Tin, with associated kernel 

pNT(e,a,b) 

a N (x,y) :=t, N 5 y (x) =T(e,a,by 1 / G(x,t;y)dt, 

JT(e,a,b) 

where G is the Green function associated with S and described in Proposition 
13.71 By Remark 15.51 an d the Green function bounds (13.351) . the continuization 
error N := T, N — lt N is uniformly bounded as an operator from X 2 — > Xi 
and, for || ■ \\x 2 uniformly bounded, is uniformly convergent from £>2 — > Bi, 
with limit uniformly Lipschitz with respect to (a,b). Thus, we may discard 
in what follows, and work directly with continuization Ejv- 
Following (13.271) . decompose 

(5.15) Sjy = + $N + "T^Ni 

with kernels E N := Sj^S y (x), Sn := Sn5 v (x), Rn '■= Ti-N&y{%) defined by 

pNT(e,a,b) 

E N '■= T(e, a, b)^ 1 / £(e, x,t;y)dt, 

JT(e,a,b) 
pNT(e,a,b) 

(5.16) Sn := T(e, a, b)^ 1 / S(e, x 1 t;y)dt, 

JT(e,a,b) 
rNT(e,a,b) 

Rn '■= T(e,a,b)' 1 / 1Z{s ) x 1 t]y)dt 1 

JT(e,a,b) 

where £, S, 71 are as in (I3.27p . 

The operators £n and Sn may be estimated as were Jn and Kn in 
Sections 15.1.11 and 15.1.21 respectively. Indeed, £n is exactly a superposition 
of terms of form j N) while Sn is a superposition of approximate Gaussian 
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terms obeying the same estimates used to bound /Cat, so can be handled in 
the same way; see [Rj IHRl IHRZj for similar calculations. 

We now turn to the residual 1Z n- By bounds f)3.32p and (13.331) . lZ y may 
be decomposed into terms 



0(e~" n ^)t~ 1 e~ ( - x ' y ~ a ^ t)2/Mt 

k=l 

of order \J y \, J as defined in Section 15.1. 2\ terms 

n 
k=l 

+ E x w >| y | } o((« + i)-^r 1 ) c -t--^m* c ^ 

a~>0, aj<0 
a">0,a+>0 

of order l-K^I, a neglible term of order e~^\ x ~y\ +t )^ and the time-derivative 
d t r(x, t; y) of terms 

r{x,t; y) =0(e- vlyl (t+l)- 1/2 ) 

^ e -(x-y-a-t) 2 /Mt + e -{x-y-a+tf/Mt\ 



of order \K\. 

Terms of order | J y \ may be estimated exactly as was J y in the proof of 
Lemma I5.6[ since those arguments depended only on modulus bounds. Lik- 
wise, terms of order \K yy \ may be handled by the modulus-bound arguments 
used in the proofs of Lemma I5TT1 and Lemma [53 to bound terms of the same 
order. Terms of order e~ v ^ x ~ y ^ +t ^ may be handled by a similar argument, us- 
ing | e~ v( -\ x ~ y \ +t ^ dt\ < C e^^ 1 - e~"^ x ~ y y Finally, derivative terms may be 
estimated by a cancellation argument like that used to bound K y , integrating 
in time to obtain fj^ r +pS)T d t r(x, t; y)dt = r(x, t; y)\ < ^ T i ~ P ' )T , hence 

r(N+ P )T 

T- 1 d t r{x,t;y)dt <C{NT)- l l\ 

J NT L2 ( x ) 
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giving convergence in C(B 2 , £>i), uniformly in T. Likewise, the limiting kernel 

r+oo 

T- 1 J d t r(x, t; y)dt = -T^rfo T;y) = {e^^ 2 ) 

is clearly bounded by C(l + \x — hence the limiting operator belongs 

to C(X 2 ,Xi) by the argument used in Lemma 15^41 to treat K y . 

Thus, the sequence TZ^ is convergent in C(X 2 , Bi), with limiting operator 
bounded in C(B 2 , B\) and C(X 2 , Xi). With Proposition 12.51 Lemmas 15.41 and 
15.61 and the above remark on the discretization error, we may conclude that 
S satisfies Assumption 12.41 

We now check that Assumption 12.241 is satisfied. The results of Lemmas 
15.41 and 15.61 clearly extend to model scattered and excited terms 

(5.17) K £ (x, t; y) := r V2 e -(*- y -a( e )t) 2 /« ? 

(5.18) J £ (x,t;y) := (w £ )'(x)errfn((-?/ - a(e)t)/2t 1/2 ), 

now including e-dependence, as long as a(e) is bounded away from 0, with 
uniform convergence with respect to e, for e in a neigborhood of the origin. 

Returning to the actual Green function, notice that e-dependence of S, 
S and 1Z is similar to (I5.17P and (I5.18p . with additional multiplicative terms 
(defined in terms of the spectral elements , lj , r±) that do not modify the 
estimates of Section I5TT1 Besides, S, S and 71 depend on a, b only through T, 
the Lipschitz function given by Lemma [4.41 Thus, convergence in C(X 2 ,Bi) 
of En + Sn+Rn is uniform with respect to e, a, b (with b G X\, \e\ + \a\ + \\b\\x 1 
small enough). Together with continuity of the Green function with respect 
to e (a consequence of bound (I3.17P ) and of S with respect to T, this implies 
continuity of the right inverse (Id — S)^ 1 with respect to e, a, b. 

Lipschitz continuity with respect to b of is straightforward, simi- 
larly as in the proof of Lemma I5.4[ since 6-derivatives of limiting operators 
T _1 pass onto the factor T~ x and the lower limit of integration only. 
Thus, Lipschitz continuity of (Id — S)^ 1 with respect to b G B\ in norm 
£(B 2 , Bi) follows from Lemmas 15.41 and I5.6[ and from Lipschitz continuity of 
T in b G B x . 

Finally, the uniform bounds (12.141) follow from the corresponding uniform 
bounds proved in Lemmas 15.41 and 15.61 and the above description of 71. □ 
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Consider now the Poincare return map system (14.281) associated with the 
truncated equations (14.171) . We can apply Proposition 12.281 for instance 
with Q = : there exists a function a —>■ e(a), defined for small \a\, such that 
e(0) = 0, e(-) continuous at 0, and periodic maps r, 9, v such that (r, 9, v) is 
a solution of (14.171) . with r(0) = a, v (0) = B(e(a), a, 0), with period 



T*(a) := T(e(a), a, B(e{a), a, 0)), 

where T is given by Lemma 14.41 and B by Proposition 12.281 From (12.171) 
and (I4.18p (ii). we find that if a is small enough, (14. 19j) is satisfied. That is, 
(r, 9, v) satisfies not only the truncated system (14.171) but also the original 
system (14.61) . and 

is a solution of the perturbation equations (14. 2p with e = e(a). At time t = 0, 

u a (-, 0) = a(<P% + 01) + B(e(a), a, 0), 
so that, with (12TT7D . 

C- l a < |K(-,0)|| Xl < Ca. 

The map u a + u e ^ is a time-periodic solution of the original, unperturbed 
equations (11.201) . satisfying (ll.23p . 

Remark 5.8. As discussed in Remark 11.81 operator (/ — S)^ 1 = ^JLqS^ 
does not preserve zero mass, fudx = 0, since J2f =0 d x includes terms of 
form K(x, NT;y)rfl±, a^>0 (resp. aj < 0) with nonzero mass but zero 
L 2 -limit as N —>■ oo. However, it does appear to preserve J l £ -udx = 0, for £ £ 
orthogonal to r^ 1 , > (resp. aj < 0), though we shall not pursue this here. 
That is, it appears to preserve the invariant subspace Kerllo complementary 
to KerL, where n is the generalized spectral projection defined in Remark 
Compare the discussion of the finite-dimensional case in Remark 12.101 



5.3 Uniqueness 

We complete the proof of Theorem 11.31 in this section, as we establish Lip- 
schitz (and later C 1 ) regularity of e(-) and uniqueness up to translation of 
solutions, using still further cancellation in the e-derivative estimates to ver- 
ify Assumption 12.61 
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5.3.1 Reduction to zero speed 



Lemma 5.9. Under assumptions (H0)-(H2), (T> £ ), let u £ , \im. z ^,± 00 u £ = u±, 
be a bounded standing-wave solution of (11.201) . and u a (x — a a t,t) another 
bounded solution, withu a (x,t) time-periodic and \im z ^± 00 u a (z) — u±. Then, 
a a = ifu £ + ^ut- 

Proof. First, observe, by standard parabolic smoothing estimates, that both 
d x u £ and d x u a converge to zero as x — > ±oo. Integrating F(u £ , e) x — u £ xx = 
over [— M, +M] for M > sufficiently large, we thus obtain the Rankine- 
Hugoniot conditions 

(RH) F{u £ + ,e) -F(u E _,e) = 0. 

On the other hand, integrating u° — a a u x + F(u a , e) x — u xx = by parts over 
the rectangle (x,t) G [-M, +M] x [0,T], M > sufficiently large and T the 
period of u a , we obtain 

(aRH) F(u E + ,e) - F(u £ _,s) = o a \u% - u e _) , 

yielding o~ a = unless u £ + = u £ _. □ 



5.3.2 Reduction to a wedge 

Recall that the analysis of Section H] applies to solutions in a wedge < 
C\a\, (a,b) := (w,v)\ t =o, with period T = T(e,a,b) determined by the 
Poincare return map construction. We must therefore eliminate the pos- 
sibility that there may exist other periodic solutions outside the wedge or 
with other periods T. 

Accordingly, we look for periodic solutions of period T of system (14. 6p . 
where T is now considered as an arbitrary parameter. The Poincare return 
map for this system is denoted by (f,g); periodic solutions of (14. 6 p are in 
one-to-one correspondence with solutions of the system 

= fie, a, b), 

(5.19) 

= g{s,a,b) = (S(e, T) - Id) 6 + N 2 (e, a, b, T), 
(a 1 i) = Ki))| 1=0 6l 2 xI ll 

(5.20) S(e,T) := e L{£)T fl £ . 
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Lemma 5.10. Under Assumptions (H0)-(H2), (T> £ ), on subspace X 2 C £>i 
containing 'Ran.ge{N 2 \x 1 ) , (Id — S(e, T)) has a right inverse (Id — S(e, T")) -1 , 
i/iai zs bounded in £(X 2 ,Xl) and £(^2, Bi), uniformly in (e,T) for T G 
[T ,Ti], < T < Ti < oo, and \e\ small enough. 

Proof. Identical to the proof of Lemma 15.71 in Section 15. 2\ which uses only 
the form of e^^IF and the fact that T is uniformly bounded from zero and 
infinity. □ 

Lemma 5.11. Under Assumptions (H0)-(H2), {T> e ), the linear operator 
S(e,T) defined in ( 15.201) satisfies A ssumption 12 . 1 51 

Proof. Consider the family of maps 

$ £ :Rx(K 2 x X x ) — ► K 2 x X h 

where $ e (c, a, b) is defined as 

x i— > (a, b(x + c)) + (u £ (x + c) — u £ (x)y 

Decomposing u e into the sum of IF-u 5 e M? and n £ M e e X\, we see that 
$ £ (c, a, 6) does belong to I 2 x li. 

Let (a, b) be an equilibrium (that is, a solution of (15.191) ) for a given 
value of e; let w e = ty + t> be the time-periodic solution of the perturbation 
equations (14. 2ft issued from a + 6, and let 

(x, t) -u e (x, t) = u £ (x + c, t) + (-u £ (a; + c) — u £ (x)) . 

The map u £ + m £ is time-periodic, and a translate of w e + u £ , itself a solution 
of the autonomous system (11.201) . Hence u £ + u £ is a solution of (11.201) , and 
u £ is a time-periodic solution of the perturbation equations (14. 2 p . The initial 
datum for u £ is $(c, a, 6). Thus, for any c, $ e (c, ■, ■) preserves equilibria. 

It is obvious that $ e is a group action on the set of equilibria. The 
derivative of $ e in its first variable at is 

9x^(0, 0,0) = (u £ )'(x), 

and it is a consequence of (T> e ) that (u 6 )' generates ker(Id — S). Finally, 
$ £ (c, a, b) is continous in c by density in X\ of C°° nXj. □ 
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Corollary 5.12. Under Assumptions (H0)-(H2), {T> £ ), any time-periodic 
solution u £ + u £ of (11.201) with ||w e ||xi sufficiently small, has a translate 
(u £ + u £ )(x) = (u £ + u £ )(x + c) such that (a,b) : = u £ (t = 0) G S e x t £ 
satisfies ||&||xi < C\a\, for any C > 0. Thus, in seeking periodic solutions, it 
is sufficient to look in the wedge \\b\\x 1 < C*i | a | . 

Proof. Just like the quadratic term A*2 in Lemma 14.51 the quadratic term 
N2 in (I5.19P satisfies (12.41) (i). Indeed, the change of coordinatization in the 
plane does not modify the Duhamel estimate. Besides, Lemma 
15.101 states that the operator S(e,T) satisfies Assumption 12.41 (up to a slight 
change in notation; indeed here S depends on the real parameters e, T, when 
Assumption 12.41 deals with operators depending on e, a, b G R x IR n x X%\ 
it is readily seen that this makes no difference in the proof of Proposition 
I2.16p . Lemma [5.111 states that Assumption 12.151 is also satisfied, and we can 
conclude by Proposition 12.161 

□ 



5.3.3 Lipschitz regularity: model computation 



Restricting to zero-speed solutions within the wedge, we now investigate the 
Lipschitz regularity of e(o). We first show how to verify Lipschitz regularity 
in e for the model problem 



(5.21] 



G £ (x,t;y) = K e (x,t;y) + J £ (x,t;y), 



now including e-dependence, where K £ is given by (15.171) and J £ by ()5.18|) . 
where a(e) < a < and a(-) is C 1 in e. By inspection, this is faithful to the 
approximation G ~ £ £ + S £ obtained by neglecting residual 1Z £ in (I3.27p . 

By Remark 1 5. 5 [ it is sufficient to treat the continuous approximants /Coo, 
j7oo- The J-term may be shown to be Lipschitz in e by a modulus-bound 
computation, 



\d e (TJ^)d x f\\ L2 < C / J £ d x f{y)dy 



L 2 



(5.22) 



+ 
+ 



(d £ d x u £ )K £ (0,t; y)f(y) dydt 
(d x u £ )d £ K £ (0,t;y)f(y)dydt 



L 2 



L 2 
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The first term in (I5.22p is bounded by C||.f the second term is bounded 
by (711(9^/11x2? due to the uniform bound IIc^c^m 6 !!^ < C, and using d e K e = 
(d e a)tKy, we see that the third term is bounded by C||9 x /||x 2 as i n fl5.11j) - 

The i^-term by (15.91) may be expressed as 
(5.23) 

d^aTt^dJ = - J d e K%x,t;y)f(y) dy + (d e T) J K £ yy (x,t;y)f(y) dy 

d e K £ yy (x,t;y)f(y) dydt. 

Lemma 15.11 implies that the L? norm of the first two terms in the right-hand 
side of (15.231) is bounded by C||/||z,°o. The third term may be bounded by 
a cancellation estimate like that used to bound in the first place, using 

d e K w = ( d e a ) tK yyyi K yyy = ~ a ~ lK yyt + ^ K yyyy and integration by parts 
to obtain 

a / K e yyy {x, t; y)t dt = K e yy {x, T; y)T + / K £ yy (x, t; y) dt 

JT(e,a,b) JT(e,a,b) 

r+oo 

+ / K yyyy(Xit;y)tdt, 
JT(e,a,b) 

where the last two terms by \K yyyy \t ~ \K yy \ are similar order and uniformly 
absolutely convergent in L 2 (x). 



5.3.4 Lipschitz regularity: full computation 

With these observations, it is straightforward to treat the full problem. 

Lemma 5.13. Under Assumptions (H0)-(H2), (T> e ), the linear operator S 
defined in (I4.29P satisfies A ssumption 12 , 61 

Proof. That S satisfies Assumptions 12.41 and 12.241 was checked in Lemma IBTTl 
It remains thus only to prove that (I2.13P (ii) is satisfied. We check below that 
^JIqS 1 - 7 is Lipschitz with respect to e in £(X 2 ,£>i). The dependence of S 
on a is through T = T(e, a, b), and so Lipschitz regularity of Yl'jLo ^ with 
respect to a is proved in the same way. 

By Remark [575] and Green function bound (I3.43p . we can neglect the error 
made in the approximation of := ^2f=i by its continuization S^v, with 



83 



associated kernel 



oo 



(5.24) a{x, y) := H6 v (x) = T(e, a, by 1 / G(x, t; y)dt, 

JT(e,a,b) 

where G is the Green function associated with S and described in Proposition 
13.71 Decomposing 

t £ = £ £ + S £ + TZ £ , 
with kernels E £ := £ £ 5 y (x), S £ := <S £ 5 y (x), R £ := TZ e S y (x) defined by 



E £ 


:=T(e, 


a^b)- 1 


poo 

/ £ £ ( 


x,t;y)dt, 








JT(e,a,b) 

POO 




S £ 


:= T(e, 


a^b)- 1 


/ S £ { 


[x,t;y)dt, 








>T{£,a,b) 




R £ 


:=T(e, 


a^b)- 1 


POO 

/ K £ 


(x,t;y)dt, 








)T(e,a,b) 





where £ £ , S £ , 1Z £ are as in (13.271) . we may estimate operators £ £ and S £ as 
we did K £ and J £ in Section 15.3.31 Thus, we need only estimate the residual 
term 7Z £ , showing uniform boundedness of ||d e 7?. E || £2 (,,,). 
By the computations of Section 15.21 we have 

R £ y = T(e, a, b)' 1 ( f {K £ y - d t r) (x, t; y)dt - r(x, T(e, a, b); y)) . 

V JT(e,a,b) J 

We have d £ T bounded, by Lemma [4.41 Likewise, d £ r(x,T;y) is bounded in 
L 2 (x) by (13.411) and boundedness of \T\. Thus, we need only show uniform 
boundedness in L 2 (x) of 

(d £ n £ y - d £ d t r)(x,t;y)dt, 



'T(e,a,b) 

which, by (13.401) . (I3.42j) may be expressed as 



POO 

/ (d £ 1¥-d £ d t r - d t f) (x, t; y)dt + lim f{x, t; y)\^ = 

JT(e,a,b) M ^°° 
poo 

/ (d £ K £ y - d £ d t r - dtf) (x, t- y)dt - f(x, T;y)\, 

JT(e,a,b) 
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where ||r||x2(a;) is evidently bounded. Observing that 

(d £ K e y - d £ d t r - d t f) (x, t; y) 

by (13.401) may be decomposed into terms Y^=i^ (e~ v ^) t~ 1 e~ ( - x ^ y ^ a ^ 
of order | J y | and terms 



t) 2 /Mt 



E° ((* + l)~ 1/2 e~ vx+ ) t- x e-^- a ^f m 

k=l 

+ E ^ a - tl > m o((t+ir^)e-^ 

a k >0, a. <0 

+ E *{\a- k t\>\ y \ } 0{{t + lY l/2 r l )e-^ -»> >-e . 

a~>0,a+>0 

of order \K yy \, we find by the calculations of Sections 15. 1.21 and 15. 1. II that 

poo 

(d £ Tl £ y - d £ d t r) (x,t;y)dt, 



-z+) 2 /Mt- V x- 



T{e,a,b) V ■ £2(3.) 

is uniformly bounded as well, completing the proof. □ 

5.3.5 Conclusion: uniqueness up to group invariance 

We now sum up the results of Sections 15.3.21 to [5.3.41 

Consider a periodic solution u e of the perturbation equations 114. 2B . By 
Corollary 15.121 u e has a translate u £ whose initial condition u £ (t = 0) = 
((a cos 6*o, a sin 6*0), b) 6 E £ x E £ satisfies 

(5.25) \\b\\ Xl <C\a\. 

If a = 0, then it 6 , hence u £ identically vanishes, by uniqueness of the solution 
of (14.21) . If a 7^ 0, coordinatize u e by (w,v) as Section HJ [w,v) solves (14. 6p . 
Consider now the associated truncated system (14.171) . where ip is defined by 
fT4TT2l) with C := C. (C as in f!Q5] .) It is a consequence of (I5T25D that (w, v) 
also satisfies system (14.171) . We can thus apply Proposition 14. 2| in particular, 
w never vanishes and we can use polar coordinates w = (re l9 , re~ l6 ), without 
loss of generality taking #(0) = 9 = 0, w(0) = (a, 0). By Lemma H~4"t 
6(T(e,a,b)) = 2tt. By Lemmas 15.131 and 15.111 the assumptions of Corollary 
(I2.17P are satisfied, and we can conclude that, up to group invariance, e = 
So(a) and b = B(eo(a), a, 0), with the notations of Proposition 12.121 



85 



5.4 C dependence 

We have so far carried out the analysis entirely in the Lipschitz framework 
of Section [21 thereby obtaining Lipschitz dependence of functions e(-), T(-), 
vS' on a. However, it is straightforward to improve this regularity to C 1 . 
Focusing on the special solution u = of Proposition 12. \7\ notice that 

= 11^(6,0,0)11* <C|a| 2 . 

Thus, we may replace linear truncation (14. lip by a quadratic order truncation 

v ■= i/j(C\r\ 2 /\v\)v 

without affecting the validity of truncated equations (14.171) . With this mod- 
ification, it is readily calculated that N r and N$ become, respectively, C 2 
and C l . In particular, we obtain from this a C l rather than Lipschitz pe- 
riod function T(e,a,b) and thereby a C l return map system (14.281) . This 
improved regularity may be carried throughout the bifurcation analysis to 
yield the result. 



6 The multi-d case: Proof of Theorem 1.4 



We now briefly describe the extension to the multidimensional case, model 
problem III, Section [T31 Consider a one-parameter family of standing planar 
viscous shock solutions u e (xi) of a smoothly- varying family of conservation 
laws 

d 

(6.1) u t = F(e,u):=A x u-J2 Fj ( £ ,u) Xj , ueW 1 

3=1 

in an infinite cylinder 

C := {x : (xi,x) G M 1 x Q}, x = (x 2 , . . . , Xj) 

Q G IR d_1 bounded, with Neumann boundary conditions 

du/dx ■ Uq = for x G <9f2, 

(or, in the case that Q is rectangular, periodic boundary conditions), with 
associated linearized operators 

d 

(6.2) L(e) :=dF/du\ u=U e = - 9 Xj A j (x 1: e) + A„ 
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A j (x,e) := Fl(u £ (x),e), denoting A 3 ±(e) := lim^ ±00 A j (z, e) = F v (u±,e). 
Profiles u 6 satisfy the standing-wave ODE 

(6.3) u' = F 1 (u,e)-F 1 (u-,e). 

Following [Z4\ |ZT] , assume: 
(HO) F j eC k ,k> 2. 

(HI) a(A±(e)) real, distinct, and nonzero, and o"(X! 0^±( £ )) rea ^ 
and semisimple for £ G M. d . 

(H2) Considered as connecting orbits of (16. 3p . u e are transverse and 
unique up to translation, with dimensions of the stable subpace S(A l + ) and 
the unstable subspace U(A X _) summing for each ston + 1. 

As in Remark 13.11 (H2) implies that u e is of standard Lax type. 
6.1 Rectangular geometry 

For clarity of exposition, we specialize now to the case of a rectangular cross- 
section Q with periodic boundary conditions, without loss of generality 

n = T d . 

This case closely resembles that of the whole space, making the analysis 
particularly transparent. In particular, we may take the discrete Fourier 
transform in transverse directions x = (x2, ■ ■ ■ , Xd) to obtain a family of 
linearized ordinary differential operators in xi, 

d 

(6.4) L^e) := L (e) ^< r W.r ; . ) - |£| 2 , 

indexed by £ G Z d-1 , where L (e) := 8%. — d Xl A 1 (xx, e) is the one-dimensional 
linearized operator of (13.31) . and £ = (£ 2 , • • • are the frequency variables 
associated with coordinates x = (x 2 , ■ ■ ■ , Xj). 

Associated to each L^(e), we may define an Evans function -D e (£, A) as 
in the one-dimensional case, whose zeroes correspond in location and multi- 
plicity with eigenvalues of L^. Moreover, the eigenvalues of L consist of the 

union of eigenvalues of L| for all £ G Z d_1 , with associated eigenfunctions 

W(x) = e^' x w(xi), where w is the eigenfunction of L|. 
To (H0)-(H2) we adjoin the Evans function condition: 
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{V £ ) On a neighborhood of Z^" 1 x {KA > 0} \ {0, 0}, the only zeroes 
of .D(£, •) are (i) a zero of multiplicity one at (£, A) = (0, 0), and (ii) a crossing 
conjugate pair of zeroes A±(e) = 7(e) +ir(e) of some L^ t , with 7(0) = 0, 
9 e 7(0) > 0, and r(0) ^ 0. 

Lemma 6.1. Conditions (H0)-(H2) and (T> £ ) are equivalent to conditions 
(P)(i)-(iii) of the introduction (Section 11.2.3)) together with F G C k , k > 2, 
simplicity and nonvanishing of a{A\_{e)), semisimplicity of a(^2£jA±(e)) 
and the Lax condition 

dim SiAKe)) + dim *7(Ai(e)) = n + 1, 

with u £ (linearly) stable for e < and unstable for e > 0. 

Proof. Essentially the same as that of Lemma 13.21 but in the final asser- 
tion (here, just a comment) substituting for the one-dimensional linearized 
stability analysis of [ZH, MaZ3] a "one-and-one-half dimensional" stability 
analysis like that used below to verify Assumption 12. 4( see Remark 16.71 □ 

Remark 6.2. Eigenvalues crossing at transverse wave- number = cor- 
respond to the one-dimensional case considered in Sections hence the 
multidimensional subsumes the one- dimensional analysis. Such crossings cor- 
respond to longitudinal "galloping" or "pulsating" instabilities. Eigenvalues 
crossing at nonzero wave number correspond in the rectangular geometry to 
transverse "cellular" instabilities as discussed in [KSJ. 

Introduce Banach spaces E\ = L 2 , B2 = d x L l H L 2 , 
X 1 = {f: \f(x)\ ^(l + H)- 1 }, 

and 

X 2 = d x {f: \f{x)\ <C7(l + |xi|)- 2 }nX 1 , 
equipped with norms ||/|| Bl = \\f\\ L 2, \\d x f\\ B2 = + ||^/|| L 2, 

= ||(1 + ImD/IlL-, and \\d x f\\ X2 = ||(1 + IxiDVllL- + R/lk, 

where d x is taken in the sense of distributions. By inspection, we have that 
£>2 C Bi, X 2 C Xi, Xi C Bi, X 2 C B 2 , and the closed unit ball in X\ is 
closed as a subset of B\. 
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6.2 Linearized estimates 



Lemma 6.3. Under Assumptions (H0)-(H2), (T> £ ), associated with eigen- 
values A±(e) of L(e) are right and left eigenfunctions <p £ ± = e l ^*'^w±(xi) and 
(f) e ± = e l ^*'^w±(xi) G C k (x,e), k > 2 as in Assumption (HO) from Section 
0, exponentially decaying in up to q derivatives as x\ — *> ±oo ; and L{e)- 
invariant projection 

j=± 

onto the total (oscillatory) eigenspace S e := Span{<^} ; bounded from L q 
or B2 to W 2,p fl X2 for any 1 < q,p < 00. Moreover, cf) e ± = d x &±,r with 
$ £ G C k+1 exponentially decaying in up to k + 1 derivatives as x — *> ±00. 

Proof. Essentially identical to that of Lemma 13.51 □ 

Defining IF, E E L(e) as in Section [3T2l denote by G as in ( 13.101) (where 
x is now multi-dimensional) the Green kernel associated with the linearized 
solution operator e Lt of the linearized evolution equations u t = L(e)u, and 
by G as in (13.111) the Green kernel associated with the transverse linearized 
solution operator e^^IF. Evidently, G = O + G, where O is defined as in 

(TO. 

Lemma 6.4 (Short-time estimates). Under Assumptions (H0)-(H2), (T> e ), 
for < t < T, any fixed T > 0, and some C = C(T), 

(6.5) H^H/lk < \\e L ^flf\\ Xl < \\f\\ x „ 

(6.6) \\e L ^%f\\ B2 , ||e Z ^n^./|U 2 < Cr l l 2 \\f\\ L ^. 

(6.7) \\e L ^d x f\\ X2 , \\e l ^Wd x f\\ X2 < Cr 1 ' 2 ^! + M) 2 /IU~, 
and 

(6-8) ||d £ (e z ^rr)/|| Bl < C||/|| Bl , 

(6-9) \\d e {e L ^W)d x fh 2 < C(\\f\\ L . + \\f\\ Bl ). 

Proof. Essentially identical to that of Lemma 13.61 □ 
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Note that the Fourier transform completely decouples the linearized prob- 
lem in the sense that the Fourier-transformed operators 

G i {xi,t;y 1 ) = FGF- 1 

acting in frequency space, T denoting Fourier transform, have simple form 

(6.10) 4 e ^> (ll) l + 



IF; := Id — IF|, where 



is the one-dimensional projection onto the oscillatory subspace of L^. 

Proposition 6.5 (Global bounds). Under Assumptions (H0)-(H2), (T> e ): (i) 
Go satisfies the pointwise bounds stated in Proposition \3. 7[ (ii) For £ 7^ 0, 
i/ie Fourier transformed solution operators e 1 ^ 1 for £ 7^ £* ; e^^^n^ /or 

£ = £* satisfy for some C, rj > 0, the exponential bounds 

11(1 + kiDe^VlU^) < <7e-*||(l + I^D/IUoc^), 
11(1 + lariDe^^/Hi-^) < Ce^||(l + N)/|U~ (:ci) . 

Proof. Assertion (i) follows immediately from the observation that Lq is ex- 
actly the one-dimensional linearized operator studied in Section [3j together 
with the fact that the bounds on the restriction to £ in case = are 
the same as the bounds on the full solution operator in the case £ 7^ of 
linearized one- dimensional stability. 

Assertion (ii) follows for |£| large by standard semigroup/asymptotic ODE 
estimates. For |£| bounded but nonzero, it follows by standard semigroup 
estimates, together with the assumption that there are no eigenvalues of 
other than the crossing pair at £ = and the computation as in fll.9p for the 
one-dimensional case that cr ess (L|) C {A : 3?A < — 2r]|£| 2 } for some rj > 0. □ 

Remark 6.6. The argument for (ii) of course fails at £ = 0, due to the lack 
of spectral gap. 
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6.3 Proof of Theorem ITU 



We restrict for definiteness to the case of periodic boundary conditions. The 
proof in the Neumann case is essentially identical. 

Using the bounds of Lemmas 16.31 and I6.4[ we may carry out the polar 
coordinate and Poincare return map construction of Section H] essentially 
unchanged to obtain a multidimensional version of Lemma 14.51 (statement 
unchanged), reducing the problem to the abstract form studied in Section 
|2j Thus, it is sufficient to establish the linearized estimates of Assumptions 
12.41 and 12.61 But, these follow easily using the one-dimensional bounds of 
Proposition I6.5( i) and the calculations of the one- dimensional case together 
with the £ ^ bounds of Proposition I6.5( ii). 

For example, uniform boundedness of J^°° e L ( e ^dt from X 2 — > X\ follows 

by 



(6.12) 



e L{£)t f dt 



T 



Xi 



l + |xi|) / e L(£)t fdt 



< 



(n 


-N) 


L 









L°°(xi,x) 

e L t {£)t fdt 



L°°(xi) 



and the observation, by Proposition I6.5( ii). Hausdorff-Young's inequality, 

and the fact that L°° controls L 1 on bounded domains, that 

(6.13) 

<T( Ce^dt 



X 



IKl + lariD/O.Qll^i) 
< C 2 sup ||(1 + |xi|)/(.,|)|| L .»( a!l) 



< C 2 ||(l + \xi\)J WLi&L^ix!)) 

<C 3 \\(1 + \x 1 \)f\\ L o o(x) 

= a 



'3H J WXi, 

together with the computation, using the one-dimensional estimates carried 
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out in Sectional and denoting / = d xl Fi + d x F, so that / = d xl Fi, of 
(6.14) 



;i + N) 



M £)t fdt 



L°°(x\) 



1 + N) / e^d^Ftdt 



L°°(xi) 



<q|(i + N) 2 F(-,o)|U~ (xi) 

< C* sup ||(1 + |x 1 |) 2 F(-,0|| LOO(a;i) 

I 

< C||(l + Ixil^FHii^oo^)) 
<C 2 ||(1 + |x 1 |) 2 F|| i o 
:=C 2 ||/|k. 



'(*) 



Other computations follow similarly. 

Remark 6.7. We point out that we have in passing set up a framework 
suitable for the linearixe stability analysis of flow in a duct, a problem of in- 
terest in its own right and somewhat different from either the one-dimensional 
case considered in [ZH, MaZ3j or the multi-dimensional case considered in 
[ZS| IZlj . It would be very interesting to try to carry out a full nonlinear 
stability analysis by this technique. 



6.4 General cross-sectional geometry 

We may treat general cross-sections Q by separation of variables, decompos- 
ing 

L(e, xi, d x ) = L q {e, xi, d xi ) + M(e, x 1 , d x ), 

where 

d 

M:=-J2^(x 1 ,e)d X] +A x , 

3=2 

and expanding the perturbation u in eigenf unctions Wj(e, x±) of M(e, x\, d x ) 
on domain Q as 



u(x 



t) = ^a fc (xi)w fc (e,xi)(x) 



fc=0 



a k =: u(k), to recover a decoupled system 
d t a k = L k a k 

:= (Lq + v k )a k + 0(\d xi w k \ + \d 2 Xi w k \)a k + 0(\d Xl w k \)d x a k , 
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where = Vk{ x i) are the eigenvalues of M associated with Wk, with both 
Vk and Wk converging exponentially in x\ to limits as x% — > ±00, that is, for 
which Lk is of the same basic form as L| in Section [6T3l Thus, we can carry 
out the entire Evans function construction of Section I6TT} with k replacing £. 

More, since Wq = 1 for Neumann boundary conditions, we have again 
that Lk = L for k = 0, justifying our notation. Likewise, augmenting P(i) 
with the assumption that Vk < —v* < for all k > 1, we recover a spectral 
gap for all Lk, k > 1. From these two observations, we obtain Proposition 16. 51 
exactly as before. Lemma EH holds also, indeed was independent of Q. Thus, 
we may carry out the entire argument of Theorem 11.41 essentially unchanged, 
provided that we can establish the generalized Hausdorff- Young inequalities 

|p|k°°(fc,a;i) < CIMIl^Oi,!, 1 ^)) 

which follow, for example, if one can establish uniform supremum bounds on 
the right and left eigenfunctions of M by asymptotic eigenvalue-eigenfunction 
estimates as k — > 00. 

Alternatively, we could perform a more general analysis allowing viscosi- 
ties with cross-derivatives, and also avoiding the need for asymptotic spectral 
analysis, by expanding u instead in eigenfunctions of (or, more generally, 
the second-order elliptic operator in derivatives of x appearing in operator 
L), and noting that for Neumann boundary conditions the first eigenfunction 
Wq is still = 1. Thus, though different wave numbers no longer completely 
decouple, we still obtain a cascaded system in which the zero wave-number 
equation is again just d t wo = L wo as in the one-dimensional case, and we 
can perform an analysis as before, but dividing only into the two blocks k = 
and k > 1, treating coupling terms in the k > 1 part as source terms driving 
the equation. 

Estimating the one- dimensional part as before in (I6.14p and for the k > 1 



93 



part substituting for fl6.12p - fl6.13p the Sobolev/Parseval estimates 



e L ^U k>l fdt 



x. 



(1 + \xi 



e Li£)t fdt 



L aD (x 1 ,x) 



< 11(1 + N) / e^U^fdt 



< 



E 

k>l 



+00 2 

(l + ixiDiKi + H)^)^^/!!^)* 



L°°(a;i) 



< 



'1 + ki|)/(- ) |)IU°=( ! riy 2 (A!) 



^CalKl + lxiD/O.DlUcc^)^) 

< C 2 ||(l + M )/ 1 1 L 2 



E 

fe>l 



(l + t-(«*-i)/2)Ce-^l*dt 



T 



<C s ||(l + |x 1 |)/|| L . 



(as) = Call J ||X 



we would then obtain the result without recourse to asymptotic eigenvalue 
estimates. We shall not pursue these issues further here. 
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